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PREFACE 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


horsepower  (550  foot-pounds/sec) 
inches 


pounds  (mass) 


square  feet 


0.3048 

745.6999 

2.54 

0.4535924 

0.09290304 


To  Obtain 
metres 
watts 

centimetres 
kilograms 
square  metres 


UNIDIRECTIONAL  SPECTRAL  WAVE  GENERATION  AND 
ANALYSIS  IN  WAVE  BASINS 

PART  I:  INTRODUCTION 


Background  and  Purpose 


1.  The  directional  spectral  wave  generator  (DSWG)  at  the  US  Army  Engi¬ 
neer  Waterways  Experiment  Station's  (WES's)  Coastal  Engineering  Research 
Center  (CERC)  provides  a  means  to  simulate  real  ocean  environments  in  a  labo¬ 
ratory  physical  model  by  incorporating  directional  spreading.  Structures 
respond  differently  to  spectral  waves  than  they  do  to  monochromatic  waves. 
Thus,  the  transformation  of  wave  energy  in  shallow  and  intermediate  water 
depths  due  to  shoaling,  refraction,  and  diffraction  on  coastal  structures  such 
as  harbors,  breakwaters,  Jetties,  mounds,  artificial  islands,  and  ship 
moorings  can  be  studied. 

2.  Unidirectional  spectral  performance  tests  were  conducted  from  Jan¬ 
uary  to  April  1987.  A  series  of  36  spectral  waves  was  generated  to  determine 
the  relationship  between  the  DSWG  control  signal  and  the  measured  waves  at 
seven  different  locations  in  the  wave  basin.  The  series  consisted  of  combina¬ 
tions  of  three  wave  periods,  two  zero-moment  wave  heights,  three  spectral 
peakedness  parameters,  and  two  peak  wave  directions.  Wave  periods  were  1, 

1.5,  and  2  sec,  wave  heights  1  and  3  in.,*  peakedness  parameters  3.3  and  7, 
and  wave  directions  0  and  -30  deg.  In  addition,  two  tests  with  peakedness 
parameters  of  15  were  run  at  a  peak  wave  direction  of  0  deg  to  investigate  the 
effect  of  spectral  peakedness.  Four  tests  were  repeated  to  verify  repeat¬ 
ability.  Finally,  four  tests  were  run  at  a  peak  direction  of  30  deg  to  note 
any  differences  in  basin  response  due  to  wave  orientation  in  the  basin. 

3.  The  purpose  of  this  report  is  to  present  the  unidirectional  spectral 
wave  test  results.  It  deals  with  the  mechanics  of  accurately  simulating,  gen¬ 
erating,  measuring,  and  analyzing  unidirectional  wave  spectra.  The  verifica¬ 
tion  of  performance  characteristics  is  a  necessary  first  step  in  our  overall 
research  program  before  site-specific  studies  can  be  undertaken. 


4.  The  report  is  divided  into  two  volumes.  Volume  I  contains  the  body 
of  the  report,  presentation  of  results  for  a  representative  case,  and  an  over¬ 
all  summary  of  results.  Volume  II  is  con.Dosed  of  appendixes  which  contain 
results  for  the  rest  of  the  cases. 

5.  A  new  method  of  generating  unidirectional  spectral  waves  has  re¬ 
cently  been  implemented.  In  Part  II,  descriptions  of  this  theory  of  simula¬ 
tion  and  generation  of  unidirectional  spectral  waves  are  presented.  The  gen¬ 
eration  program  simulates  a  stroke  time  series  for  each  of  the  6l  independent 
paddles  of  the  DSWG  using  frequency  domain  methods  and  a  "235’’  fast  fourier 
transform  (FFT)  algorithm.  Either  deterministic  amplitude  and  random  phase 
(DSA)  or  random  amplitude  and  phase  (NSA)  options  are  available.  Two  cases 
were  run  to  compare  these  two  simulation  methods.  The  unidirectional  wave 
spectrum  is  simulated  as  a  product  of  a  frequency  spectrum  and  a  directional 
spreading  function.  The  desired  frequency  spectra  may  be  an  empirical  TMA  or 
Ochi-Hubble  or  discrete  input  at  any  selected  frequencies.  The  directional 
spreading  function  may  be  an  empirically  derived  wrapped  normal  or  discrete 
input  at  equally  spaced  angles  in  a  360-deg  circle.  Although  the  directional 
spectrum  muy  to  mult  i.'i — ir.  both  f.cquer.wy  cr.d  ui.cction  space,  only  uni- 
modal  directional  spectra  were  tested  in  this  researen  effort.  Also,  only  the 
TMA  frequency  spectrum  and  the  wrapped  normal  spreading  function  were  tested. 

6.  Part  II  also  contains  a  description  of  the  linear  transfer  function 
for  directional  "piston  mode"  wave  makers  ror  conv  *  Cing  a  «u*t  c  ltva’u  iuu  kj  j.iue 
series  to  corresponding  strokes.  These  values  are  converted  to  digital 
voltage  units  prior  to  signal  generation. 

7.  Part  III  describes  the  theory  of  time  and  frequency  domain  analyses 
of  the  measured  data.  Included  are  zero-downcrossing,  spectral  and  unidirec¬ 
tional  spectral,  frequency  response  functions,  and  wave  grouping.  Three  dif¬ 
ferent  methods  of  analysis  are  available  for  estimating  the  unidirectional 
spectral  density  estimates:  (a)  an  eigenfunction  procedure,  (b)  a  method 
based  on  the  Fourier  series  expansion  of  the  directional  spectrum,  and 

(c)  Maximum  Liklihood  Method  (MLM) .  Only  the  second  method  is  addressed  in 
this  report.  The  response  amplitude  operator  (RAO)  transfer  functions  are 
time  domain  equivalents  of  the  ratios  of  the  measured  to  predicted  fre^  '-.cv 


means  of  adjusting  the  control  signals  for  more  faithful  reproduction  of  the 
desired  spectral  shapes.  Finally,  comparisons  between  Goda's  and  Kimura's 
methods  are  used  to  provide  wave  group  statistics  for  a  run  of  high  waves  from 
the  measured  data.  Successive  wave  heights  are  assumed  to  be  dependent,  a 
fact  which  is  observed  in  nature.  Probability  distributions  and  run  lengths 
for  both  methods  are  described. 

8.  In  Part  IV,  descriptions  of  the  test  setup,  including  the  direc¬ 
tional  wave  basin  and  generator,  wave  elevation  measurement  system,  test  case 
description,  and  test  case  generation  and  analysis  parameters  are  given.  <">ne 
of  the  unique  features  of  the  DSWG  is  its  flexible-plate  seals  between  indivi¬ 
dual  paddles  which  generate  a  more  uniform  wave  form  with  less  spurious  waves. 
The  directional  wave  basin  is  lined  with  wave  absorber  frames  to  minimize 
reflections.  Resistance  or  capacitance  wave  gages  with  remote-controlled 
stepper  motors  were  used.  They  were  arranged  in  a  2-3- 1-7  linear  array  incor¬ 
porated  in  a  cross  pattern  centrally  located  20  ft  in  front  of  the  DSWG  about 
its  center  line.  The  unit  lag  spacing  was  designed  to  optimize  resolution 
within  the  lower  and  upper  cutoff  frequencies  while  minimizing  spatial 
aliasing . 

9.  Test  results  and  analyses  are  presented  in  Part  V.  Time  domain 

analysis  results  include  time  series  traces  of  wave  elevation  and  wave  height 
for  each  test  case.  Frequency  domain  results  include  comparisons  of  predicted 
and  measured  peak  wave  periods,  zero-moment  wave  height,  peak  wave  direction, 
frequency  spectra,  and  directional  spreading  function.  A  comparison  of  the 
DSA  and  NSA  methods  of  unidirectional  spectral  geueiation  is  Also 

transfer  functions  as  a  function  of  frequency  for  each  test  case  are  pre¬ 
sented.  Wave  group  statistics  for  a  run  of  high  waves  based  on  Goda's  and 
Kimura's  methods  are  presented  next.  Included  are  (a)  the  effect  of  lag 
length  on  the  correlation  coefficient,  (b)  probability  distributions  for  run 
of  high  waves  exceeding  thresholds  of  the  mean  and  significant  wave  heights, 
(c)  effect  of  correlation  coefficient  on  the  mean  run  length  for  thresholds  of 
mean  and  significant  wave  heights,  (d)  effect  of  correlation  coefficient  on 
the  total  run  length,  and  (e)  effect  of  Goda's  peakedness  parameter  on  the 
mean  run  length  and  the  correlation  coefficient. 

10.  Finally,  Part  VI  contains  a  summary  of  results  and  recommenda- 


wherever  possible  throughout  this  report  and  the  computer  software  implemented 
for  wave  generation  and  analysis. 


PART  II:  SPECTRAL  SIMULATION  AND  GENERATION  THEORY 


Directional  Wave  Spectrum  Parameterization 

11.  A  directional  or  unidirectional  wave  spectrum  S(f,0)*  is  usually 
parameterized  as  the  product  of  two  parts:  a  one-dimensional  frequency 
spectrum  S(f)  and  a  directional  spreading  function  D(f,9)  as 

s(f,e)  =  s(f)  D(f,e)  (D 


where 

f  =  frequency 
0  =  wave  direction 

The  frequency  spectrum  and  spreading  function  are  subject  to  the  following 
constraints : 


S(f)  =  / 
0 


2tt 


S(f,0)  dO 


/  D(f,e)  do  =  1  (2) 

0 

Frequency  spectrum 

12.  The  desired  or  target  frequency  spectrum  may  be  input  in  one  of 
three  ways.  It  can  be  calculated  according  to  the  formula  for  an  empirical 
TMA  or  Ochi-Hubble  or  input  discretely  from  a  numerical  or  field  observation. 
Any  number  of  frequencies  can  be  used,  and  the  spectrum  can  be  multimodal 
(i.e.  multiple  peaks).  For  the  research  reported  here,  a  unimodal,  depth- 
limited  TMA  spectral  form  was  used  since  it  is  :,iost  representative  of  inter¬ 
mediate  to  shallow-water  wave  conditions. 

13.  The  TMA  spectrum  is  a  function  of  five  parameters:  peak  frequency 
fp  ,  an  alpha  constant  a  ,  peak  enhancement  factor  y  ,  spectral  width 
parameter  o  ,  and  the  water  depth  h  .  The  first  four  parameters  are  those 


associated  with  the  Joint  North  Sea  Wave  Project  (JONSWAP)  spectrum.  Vincent 
( 198*0  showed  that  the  value  for  a  can  be  estimated  from  the  wave  steepness 
e  such  that  it  is  valid  for  all  water  depths  as  follows: 

a  =  I6w2  e2  (3) 


where 

e  =  HmO^4  V 

=  zero-moment  wave  height  = 

mg  =  area  under  the  spectrum  within  lower  and  upper  cutoff  frequencies 
f  and  f  ,  respectively 

L  =  wavelength  associated  with  fp 

The  peak  enhancement  factor  controls  the  peakedness  of  the  spectrum  and  typi¬ 
cally  varies  between  1  to  3-3  for  sea  conditions  to  7  and  higher  for  swell 
waves.  Both  a  and  y  influence  the  energy  contained  in  the  wave  spectrum. 
Typical  values  for  the  left  and  right  ob  spectral  width  parameters  are 

0.07  and  0.09,  respectively.  Like  the  JONSWAP  spectrum,  a  wide  range  of 
single  peaked  spectra  may  be  simulated  by  varying  these  five  parameters.  An 
expression  for  the  TMA  spectrum  S^.^fl(f)  in  terms  of  the  JONSWAP  spectrum 
Sj(f)  is 

STMA(f)  =  SJ(f)  *  *(f’h)  (4) 

Kitaigorodski i  et  al.  (1975)  provides  a  full  expression  for  4>  which  can  be 
approximated  by 

for  fl.  <  1 

h 

for  1  <  <  2 

for  Ph  >  2  (5) 


360  deg.  The  wrapped  normal  spreading  function  was  used  in  this  study,  .'t  is 
a  function  of  the  principal  direction  and  the  spreading  standard  devia¬ 

tion  .  Both  are  linear  functions  of  frequency  f  ,  consisting  of  a  con¬ 
stant  and  slope  component.  The  Fourier  series  representation  for  the  wrapped 
normal  spreading  function  is 

L  /.  >2 

1  1  T — '  ^ 

d(f,q)  =  5-  +  -  y  exp  ■  — 5 —  °os  ■  o  (6) 

it  n  i— '  2  m 

i=1 


where 


a 

m 


°0  +  °1 
°0  +  fl1 


(f  -  fp) 
(f  -  fp) 


The  number  of  harmonic  terms  L  is  arbitrarily  selected  to  adequately  repre¬ 
sent  the  Fourier  series  of  the  directional  spreading  function.  It  can  be  of 
the  order  of  5  or  as  large  as  6.0/cj^  .  In  this  study,  both  and 

were  set  to  zero.  For  unidirectional  spectral  waves,  a  value  of  oQ  =  1  is 
used. 

15.  One  of  the  advantages  of  the  wrapped  normal  spreading  function  is 
its  high  resolution.  For  a  unidirectional  spectrum,  it  has  a  large  value  at 
,  representative  of  a  Dirac  delta  function,  and  falls  off  rapidly  for  dif¬ 
ferences  in  angle,  0  -  0^  ,  greater  than  0  deg.  After  a  7-deg  difference, 
the  value  of  D(f,o)  is  zero.  Thus,  care  must  be  exercised  in  selecting  the 
number  of  directional  increments  for  simulating  a  unidirectional  spreading 
function.  For  instance,  for  =  15  deg  and  a  directional  increment  of 

22.5  deg,  an  erroneous  value  of  D(f,o)  =  0  at  0  would  be  calculated. 

m 

Therefore,  a  minimum  of  36  increments  at  10-deg  intervals  is  necessary  to 
ensure  that  the  angular  difference  never  exceeds  5  deg.  In  this  study, 

72  increments  at  5  deg  were  used. 


Wave  Elevation  Realizations 

16.  The  psuedo-integral  model  (Borgman  1984,  Briggs  1987c)  for  the 
unidirectional  spectral  surface  elevation  time  series  n  at  paddle  location 
(x,y)  and  time  t  is 


(7) 


n(x,y,t)  =  2  f 
0 


A(f)  exp  (-i«)  exp  ( i<*> ) 


where 

A(f )  =  amplitude  function  described  below 

9  =  independent  random  phase,  uniformly  distributed  on  (0,2rr) 

<t>  =  kx  cos  0  +  ky  sin  0  -  2uft 
k  =  wave  number 

9  =  direction  toward  which  waves  travel,  clockwise  from  x-axis 


17.  The  model  described  by  Equation  7  is  a  double  summation  model 
(Borgman  1969)  since  both  frequency  and  direction  are  independent  variables. 
Although  this  method  is  the  preferred  model  to  represent  a  directional  sea,  it 
is  phase  locked  or  phase  dependent  since  different  directions  are  summed  at 
the  same  frequency.  This  phase  locking  gives  a  fluctuating  RMS  (i.e.  root 
mean  square)  which  tends  to  go  to  zero,  however,  as  the  number  of  components 
becomes  large  (Pinkster  1984).  Thus,  a  long  time  series  should  be  used.  An 
alternate  model  to  the  double  summation  model  is  the  single  summation  model. 
Tais  model  is  sometimes  preferred  because  less  computer  time  and  space  are 
required  and  possible  phase  locking  is  not  a  concern.  In  this  model,  a  single 
direction  is  randomly  selected  for  each  frequency  component.  A  procedure  for 
selecting  this  single  direction  (Sand  and  Mynett  1987)  is  to  assign  a  cumula¬ 
tive  directional  distribution  for  each  frequency  and  select  a  value  using  a 
uniform  random  distribution  U(0,1)  . 

18.  Realizations  of  the  desired  time  series  described  in  Equation  7  for 
the  specified  unidirectional  wave  spectrum  are  simulated  in  the  frequency 
domain  using  one  of  two  methods  and  then  fourier  transformed  to  the  time 
domain.  The  first  method  is  a  deterministic  amplitude,  random  phase  only 
model  (DSA)  and  the  second  method  is  a  nondeterministic  (Rayleigh)  amplitude, 
random  phase  model  (NSA)  (Tucker,  Challenor,  and  Carter  1984;  Hudspeth,  Nath, 
and  Sollitt  1983;  Isaacson  1985;  Elgar,  Guza,  and  Seymour  1984).  Both  methods 
also  are  known  by  other  names  including  the  Random  Pnase  Method  for  the  DSA 
method  and  the  Random  Coefficient  Method  for  the  NSA  method  (Sand  and  Mynett 
1987). 

19.  In  the  DSA  method,  the  fourier  coefficients  are  calculated  from  the 
target  spectrum  with  the  deterministic  amplitude  A  constrained  to  be 


A(f )  =  V t2  S(f,e)  df  de] 


(8) 


They  are  then  coupled  with  random  phases  and  inverse  fourier  transformed  using 
a  "235"  FFT  to  give  the  surface  elevation  time  series.  The  length  of  the  time 
series  N  is  the  product  of  the  integers  2,  3,  and  5  raised  to  integer 
powers,  or 


N  =  2K  3L  5M  (9) 

In  this  method,  the  simulated  spectrum  always  matches  the  target  spectrum. 

20.  In  the  NSA  method,  Equation  7  is  rewritten  in  an  equivalent  complex 
fourier  series  form.  The  complex  amplitude  A  is  defined  by 

A(f)  =  U(f)  -  i  V(f)  (10) 


where 

U(f)  =  real  amplitude  component 
V(f)  =  imaginary  amplitude  component 
Both  components  are  independent  Rayleigh  random  variables  with  zero  mean 
and  variance  S(f,e)  df  de  .  These  U  and  V  components  are  obtained  by 
first  generating  Gaussian  distributed,  zero  mean,  unit  variance  random 
variables  which  are  then  multiplied  by  the  desired  standard  deviation  (i.e. 
V[S(f,e)  df  de) ]  .  The  inverse  method  of  generating  new  random  variables  from 
a  normal  distribution  is  used  to  obtain  the  U  and  V  components.  If  the 
normal  or  Gaussian  probability  function  is  defined 

1  r X 

P<x)  =  u  =  s  / 

— <x» 

Then  a  new  random  variable  x  given  by 


can  be  generated  using  the  approximation  from  Abramowitz  and  Stegun  (1970) 


c^  +  c2t 


(12) 


where 


t= 

C0  =  2.515517  d,  =  1.432788 

c1  =  0.802853  d2  =  0.189269 

c2  =  0.010328  d3  =  0.001308 

Again,  an  inverse  FFT  gives  the  simulated  wave  elevation  time  series.  The  NSA 
method  differs  from  the  DSA  method  in  that  the  simulated  spectral  estimates  do 
not  exactly  match  the  target  spectrum  but  vary  statistically  about  the  true  or 
desired  spectral  values.  The  variance  of  the  variance  is  said  to  be  more 
realistic  than  with  the  DSA  method. 

21.  If  there  are  many  frequency  and  directional  components  in  the  wave 
train  (i.e.  of  the  order  of  1,000  components),  there  is  little  difference  be¬ 
tween  the  DSA  and  NSA  options  for  the  complex  wave  amplitude.  Both  are  asymp¬ 
totically  equivalent  by  the  central  limit  theorem  and  behave  very  nearly  as  a 
Gaussian  process  in  that  case  (Rice  1944,  1945;  Elgar,  Guza,  and  Seymour 
1985).  However,  with  some  increase  in  computer  time,  the  addition  of  the 
Rayleigh  variable  in  the  complex  amplitude  produces  an  exact  Gaussian  process. 
Addition  of  the  Rayleigh  variable  is  particularly  useful  for  wave  trains  which 
are  very  narrow  banded  in  both  frequency  and  direction. 

22.  In  summary,  the  least  random  (i.e.  most  deterministic)  simulations 
involve  the  single  summation  model  and  the  DSA  method.  The  double  summation 
model  and  NSA  methods  produce  greater  variance  in  the  variance  of  the  synthe¬ 
sized  record.  The  procedure  used  in  this  study  is  predominantly  the  double 
summation  model  and  the  DSA  method  of  frequency  domain  simulation. 

Transfer  Functions 

23.  Once  the  wave  elevation  time  series  have  been  simulated  for  each  of 
the  61  DSWG  paddles,  they  are  converted  to  corresponding  stroke  time  series 
using  a  height-to-stroke  transfer  function.  The  three-dimensional  form  F^  , 
which  includes  directional  effects  (Sand  1979),  was  originally  derived  by 
Biesel  (1954)  and  is  valid  for  all  nondimensional  water  depths  kh  and  is 
expressed  as 


4  \ 


—  —  l-  /  <M.U 


Figure  1  illustrates  how  varies  for  frequencies  between  0.1  to  2  Hz  for 

the  two  wave  directions  of  0  and  -30  deg  used  in  this  study.  The  wave  break¬ 
ing  limit  at  the  higher  frequencies  is  also  shown.  This  highest  wave  which 
can  be  generated  at  the  wave  maker  (Ahrens  1986*)  is 

=  0.1  L  tanh  kh  (15) 


where  L  is  wavelength. 


WAVE  FREQUENCY,  HZ 

Figure  1.  Height-to-stroke  transfer  function,  including 
wave  breaking  limits 

24.  Finally,  the  stroke  time-history  is  converted  to  a  voltage  time- 
history  using  the  conversion  factor  for  the  wave  maker.  For  the  maximum 
stroke  of  12  in.,  the  DSWG  has  a  peak-to-peak  voltage  of  20  v  in  4,096  digital 
units.  Thus,  the  conversion  factor  is  341.33  A/D  units  per  inch  or  a 
resolution  of  0.00293  in. 


*  Personal Communication,  January  1986,  John  P.  Ahrens,  Research  Oceanog- 
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25.  In  all  analyses,  it  is  assumed  that  the  wave  elevation  time  series 

is  a  discrete,  real-valued  time  sequence  with  equal  time  intervals  t  =  nAt  . 

The  relationship  between  the  total  length  of  the  time  series  Tp  and  the 

total  number  of  points  N  is  Tp  =  NAt  .  The  value  of  N  ,  although  even, 

does  not  have  to  be  a  power  of  2.  The  basic  frequency  increment  is  defined  by 

Af  =  1.0/T  . 

r 


Time  Series  Analysis 

26.  The  time  series  analysis  (TSA)  package  performs  both  time  and  fre¬ 
quency  domain  analyses  of  measured  surface  elevation  data.  Currently,  the 
package  allows  12  different  processes  to  be  run  on  the  data.  Time  domain  pro¬ 
cesses  include  time  series  strip  charts  of  the  raw  data,  zero  uperossing  or 
downcrossing,  crest  and  trough  height,  and  autocorrelation  and  crosscorrela¬ 
tion.  Frequency  domain  processes  include  autospectra  or  cross-spectra,  fre¬ 
quency  response  between  channels,  coherence,  and  Goda  reflection.  Only  three 
of  the  processes  were  used  in  this  study:  time  series  strip  charts,  zero 
downcrossing,  and  single  channel  spectral  analysis.  It  is  beyond  the  scope  of 
this  report  to  give  a  detailed  description  of  the  theoretical  development  of 
each  process  in  the  TSA.  The  interested  reader  is  referred  to  the  report  by 
Long  (in  preparation)  for  more  detailed  descriptions  of  the  algorithms. 

27.  The  strip  chart  option  plots  the  raw  time  series  data  for  one  or 
more  of  the  available  gages.  The  time  series  have  not  been  detrended.  The 
plots  are  cross-scaled  to  facilitate  readability. 

28.  The  zero  downcrossing  analysis  calculates  statistics  of  wave  eleva¬ 
tion,  wave  height,  and  period  for  the  datum  selected.  This  datum  can  be 
either  the  mean  or  an  externally  imposed  value.  For  the  surface  elevation, 
the  mean,  root  mean  square,  standard  deviation,  minimum,  and  maximum  values 
are  determined.  Wave  heights  and  periods  corresponding  to  the  difference  be¬ 
tween  the  minimum  and  maximum  values  between  consecutive  datum  downcrossings 
are  calculated.  From  these  two  arrays,  mean  and  significant  values  are 
computed.  Also,  the  total  number  of  waves  (i.e.  zerocrossings)  and  the  maxi¬ 
mum  wave  height  are  given.  Wave  heights  for  each  gage  are  saved  for  later 
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distributions  and  fitted  Weibull  distributions  can  be  determined  and  plotted. 

29.  Spectral  density  estimates  are  calculated  for  individual  gages 
after  preprocessing  to  detrend  and  window  the  time  series.  Detrend  options 
include  removing  the  mean  or  a  linear  or  second-order  trend.  Window  options 
include  10  or  50  percent  cosine  bell  or  cubic  polynomial.  The  data  are 
fourier  transformed,  band  averaged  between  lower  and  upper  cutoff  frequencies, 
and  plotted.  Measured  spectral  estimates  for  each  gage  are  saved  for  later 
use  in  calculating  frequency  response  estimates.  Three  different  methods  are 
used  to  calculate  the  peak  frequency  and  corresponding  period.  The  single 
line  method  gives  the  frequency  having  the  largest  energy  content.  The  Delft 
method  gives  the  mean  frequency  of  ail  frequencies  having  energy  content 
greater  than  or  equal  to  80  percent  of  the  maximum.  The  CERC  method  gives  the 
center  frequency  of  the  band  of  11  consecutive  frequencies  having  the  largest 
total  energy  density.  Zeroth,  first,  and  second  moments  are  estimated. 

Goda's  (1970)  peakedness  parameter  Qp  is  also  calculated.  Finally,  the 
average  period  from  the  zeroth  and  second  moments  is  estimated. 

Unidirectional  Spectral  Analysis 

30.  Three  different  methods  for  estimating  S(f,0)  from  an  array  of 
N  surface  wave  gages  are  available.  These  are  (a)  an  eigenfunction  proce¬ 
dure,  (b)  a  method  based  on  the  fourier  series  expansion  of  the  directional 
spectrum,  and  (c)  the  MLM.  Because  of  its  simplicity,  only  the  fourier  method 
is  used  here. 

Spectral  formulation 

31.  The  fourier  series  method  is  based  on  the  relation  (Borgman  1979, 
Briggs,  Borgman,  and  Outlaw  1987)  that  the  autospectra  S^Cf)  and  cross¬ 
spectra  Sjj(f)  between  all  pairs  of  wave  elevation  time  series  can  be 
expressed  as  a  linear  combination  of  the  directional  components  of  S(f,0)  at 
that  frequency.  If  the  wave  elevation  time  series  n  is  defined  as  in  Equa¬ 
tion  7,  then  a  system  of  equations  for  the  spectral  matrix  of  autospectra  and 
cross-spectra  in  terms  of  S(f,9)  is 

2n 

su(f>  =  /  S(f,0)  do 


(16) 


(17) 


(f) 


S(  f  ,0)  exp  ( j )  d0 


where 

B^j  =  kX^j  cos  0  -t-  kY^j  sin  0 

XiJ  =  Xi  '  XJ 

Yu  =  yi  -  yj 

x  =  x-axis  gage  coordinates  at  location  i  or  J 
y  =  y-axis  gage  coordinates  at  location  i  or  J 

The  cross-spectra  are  composed  of  real  co-spectra  Cjj  and  imaginary  quadra¬ 
ture  Q^j  components. 

Parameterization  of  directional  spectra 

32.  The  first  step  toward  solving  the  set  of  equations  above  is  to 
parameterize  S(f,0)  as  the  product  of  a  frequency  spectrum  and  a  directional 
spreading  function  (Equation  1).  The  autospectral  density  Sii  (f)  for  each 
of  N  gages  is  estimated,  and  then  a  combined  best  estimate  S(f)  is 
obtained  using  a  harmonic  mean,  as  follows: 


S(f) 


N 

n  S..(f) 
i  =  1 


1/M 


(18) 


The  directional  spreading  function  is  initially  approximated  by  a  truncated 
fourier  series  expansion  of  L  =  5  harmonics. 

L 

D(f,0)  =  a(f)  cos  (10)  +  b  (f)  sin  (10)  (19) 

<=1t  '  j  t  l 

1=1 


where 

aa(f)  =  real  fourier  coefficient  of  the  spreading  function 
bfc(f)  =  imaginary  fourier  coefficient  of  the  spreading  function 

Spectral  estimation  of  measured  data 

33.  The  next  step  is  to  calculate  measured  autospectral  and  cross- 
spectral  density  estimates  for  each  gage  and  each  gage  pair,  respectively.  A 


10  percent  cosine  bell  window  is  applied  in  the  time  domain,  and  the  cross- 
spectral  matrix  is  calculated  using  the  "235"  FFT. 

34.  A  Gaussian  smoothing  function  is  then  used  to  smooth  the  esti¬ 
mates.  This  procedure  is  like  "band  averaging"  since  raw  spectral  estimates 
are  smoothed  in  the  frequency  domain.  However,  it  tends  to  give  a  smoother 
transition  since  it  is  more  of  a  weighted  moving  average  in  that  overlapping 
is  used.  The  Gaussian  smoothed  line  spectra  Sffl  for  each  frequency  mAf  is 
defined  as 


S 

m 


i?-j  “J  S°-J 

J 


J 

£  w 

J=-J 


(20) 


where  Wj  is  a  weighting  function  defined  below,  and  Sm_j  is  the  raw  auto- 
spectral  or  cross-spectral  estimate  at  frequency  (m-j)  Af  .  The  weights  are 
defined  over  an  integer  number  of  spectral  lines  J  equivalent  to  ±3a 
standard  deviations  of  a  Gaussian  curve,  or 


Wj  =  exp 


(.1At)2 

2o 


(21) 


Since  the  area  under  a  Gaussian  curve  =  \(2tt)  a  ,  an  equivalent  rectangle 
having  the  same  area  has  height  of  1.0  (at  x  =  0)  and  a  bandwidth  b  defined 
by 


I  —  V  c-t-  ! 

V(2*> 

35.  For  Gaussian  smoothing  the  number  of  bands  averaged  in  each 
smoothed  spectral  estimate  is  considerably  greater  than  a  comparable  frequency 
domain  band  averaging  procedure.  In  band  averaging,  the  resolution  bandwidth 
b  is  the  product  of  the  number  of  bands  averaged  and  the  basic  frequency 
increment  Af  of  the  line  spectra.  In  Gaussian  smoothing  the  number  of 
bands  M  in  the  smoothed  average  is  given  by 


(23) 


p 


Solution  procedure 

36.  The  for  each  gage  are  inserted  into  the  left-hand  side  of 

Equation  16.  They  are  also  substituted  into  Equation  18  to  calculate  S  . 

The  Sjj  are  calculated  for  each  pair  of  gages  and  substituted  into  the  left- 
hand  side  of  Equation  17.  The  estimate  for  S  and  the  parameterized  direc¬ 
tional  spreading  function  D(f,Q)  are  substituted  into  the  right-hand  side  of 
Equation  17  for  S(f,0)  .  Thus,  for  N  gages,  a  set  of  N2  simultaneous 
linear  equations  (i.e.  N  autospectral  equations  of  the  form  of  Equation  16 
and  N(N-1)/2  pairs  of  cross-spectral  equations  of  the  form  of  Equation  17) 
can  be  solved  for  the  fourier  coefficients  of  the  spreading  function.  How¬ 
ever,  since  some  of  these  equations  may  be  theoretically  zero,  the  number  of 

p 

available  equations  is  usually  less  than  N  .  The  matrix  form  of  this  set  of 
equations  is 

{SPECTRA}  =  [TRANSFER]  *  {DIR}  (24) 

The  vector  {SPECTRA}  consists  of  the  measured  Si 1  estimates.  The  vector 
{DIR}  is  composed  of  the  fourier  coefficients  for  D(f,i)  at  frequency  f 
and  is  to  be  estimated.  The  matrix  [TRANSFER]  is  therefore 


2n 

"cos' 

"cos' 

(TRANSFER]  =  / 

0 

sin 

(Bij) 

_sin_ 

(25) 


37.  For  an  array  of  surface  gages,  Equation  25  could  be  expressed  in 
the  form  of  Bessel  functions  of  order  1  (Borgman  1969,  Borgman  and  Panicker 
1970).  Since  the  number  of  available  equations  (i.e.  N2)  is  greater  than  the 

number  of  unknown  fourier  coefficients  (i.e.  10  for  5  harmonics),  a  least 
squares  fourier  transform  method  for  numerical  integration  is  used.  The  first 
step  is  inversion  of  the  (TRANSFER]  matrix  to  solve  for  the  and 

spreading  coefficients  contained  in  the  vector  {DIR}  .  The  technique  for 
doing  this  is  based  on  the  vector  linear  regression  model: 


m  =  [X]  *  {B}  ♦  {e} 


(26) 


where 


{ Y}  =  n-component  vector  of  dependent  variables 
[X]  =  n  x  k  matrix  of  independent  variables 
{B}  =  k-component  vector  or  regression  coefficients 
{ e }  =  n-component  vector  of  prediction  errors 

The  least  square  solution  for  {B}  is  given  by  the  solution  to  the  "normal 
equations" 

([X]'  •  [X])  *  {B}  =  [X]'  «  {Y}  (27) 

where  the  prime  denotes  the  matrix  transpose.  Principal  components  regression 
is  used  to  invert  Equation  24  based  on  the  development  of  the  Moore-Penrose 
generalized  inverse  from  the  eigenvector  analysis  of  (X1  *  X)  (Searle  1982). 

38.  The  advantage  of  the  Fourier  method  is  its  simplicity.  The  disad¬ 
vantage  is  arbitrary  truncation  of  the  Fourier  series  at  L  (=5)  harmonics 
where  this  number  of  coefficients  (i.e.  2L  +  1)  is  considerably  less  than  the 
number  of  equations  available.  Stepwise  regression  tends  to  minimize  this 
truncation  effect.  Jennrich  (1977)  specifies  this  least  squares  procedure 
which  sequentially  selects  the  "best"  first  harmonic  coefficients.  Next, 
"best"  second  harmonic  coefficients  are  added,  and  another  least  squares 
iteration  is  performed  with  all  coefficients.  This  procedure  continues  with 
"best"  higher  order  harmonic  coefficients  being  added  one  at  a  time  until 
little  improvement  is  realized  after  a  least  squares  iteration. 

39.  The  use  of  eigenvalue  analysis  helps  to  overcome  the  disastrous 
effects  of  an  almost  singular  X'X  matrix  in  Equation  27.  It  provides  a 
number  of  eigenvalues  which  are  significantly  nonzero  as  an  upper  limit  on  the 
number  of  coefficients  which  can  be  estimated  with  the  stepwise  regression 
technique. 

40.  The  D(f,0)  given  in  Equation  19  results  in  a  distorted  estimate 
of  the  true  spreading  function  due  to  the  truncation.  This  estimate  is 
improved  by  fitting  a  parameterized  wrapped  normal  formula  (Equation  6)  to  the 
spreading  function  (Borgman  1979).  Assuming  that  the  lower  order  fourier 
coefficients  are  the  least  distorted  by  the  truncation,  they  are  used  in  the 
parameterized  model  to  obtain  estimates  of  the  higher  order  coefficients.  A 
new  set  of  fourier  coefficients  (both  lower  and  higher  order)  is  calculated 


using  an  iterative  procedure  until  convergence  is  reached. 


Principal  direction 
and  directional  spread 

41.  If  the  wrapped  normal  spreading  function  is  defined  as  in  Equa¬ 
tion  6,  the  full-circle  fourier  coefficients  a„  and  b  are 
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The  corresponding  coefficients  for  the  first  harmonic  a^  and  are  then 
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Therefore,  the  principal  direction  0^  at  frequency  f  is 


9  (f) 
m 


arctan 


b1  (f ) 


(30) 


The  mean  of  all  e  equals  0  . 
m 

42.  Taking  the  magnitude  of  Equation  29  and  rearranging  gives  the 
directional  spread  or  standard  deviation  c?m  at  frequency  f  ,  as  follows: 


o(f)  = 
m 


-2  In 


a1  (f)  +  b 


(f) 


(3D 


43.  At  each  frequency,  (N-1)  independent  estimates  of  each  fourier 

coefficient  are  obtained.  Various  methods,  including  a  weighted  average  or 

least  squares,  could  be  used  to  select  a  "best"  choice.  A  simple  average  is 

used.  A  similar  set  of  equations  could  be  obtained  for  0  and  o  if  the 

m  m 

second  or  higher  order  harmonic  coefficients  were  used. 


Response  Amplitude  Operator  Function 

44.  The  height-to-stroke  transfer  function  discussed  in  Part  II  is 


Because  of  leakage  around  and  under  a  "wetback"  wave  maker,  electronic  and 
mechanical  losses  of  the  DSWG,  use  of  linear  wave  theory,  and  basin  response 
characteristics,  the  desired  wave  spectrum  is  usually  not  faithfully  repro¬ 
duced.  Although  the  differences  in  spectral  shape,  peak  period,  and  wave 
height  may  be  slight,  it  is  necessary  to  correct  the  control  signal  for  accu¬ 
rate  reproduction  in  future  runs.  This  correction  can  be  done  once  or  itera¬ 
tively  over  several  runs  until  the  improvement  between  corrections  is  slight. 
Thus,  an  RAO  transfer  function  is  calculated.  Although  the  RAO  was  calculated 
for  each  case,  it  was  not  used  in  this  study  to  correct  the  control  signals. 

45.  The  raw  measured  for  each  gage  and  predicted  Sp  spectral 

estimates  are  first  Gaussian  smoothed  (see  paragraph  34)  to  a  desired  band¬ 
width  which  must  be  equal  to  or  greater  than  Af  .  The  RAO  is  then  calculated 
in  the  frequency  domain  at  each  smoothed  frequency  as 


RAOi(f)  = 


Sii<f> 

Vf) 


1,2,3, ...N 


(32) 


An  average  of  all  points  within  each  band  or  the  midpoint  in  each  band  may  be 
selected  for  each  band.  In  addition  to  RAO's  for  individual  gages,  an  average 
of  all  or  a  selected  number  of  gages  may  be  calculated.  The  RAO  is  set  equal 
to  1.0  for  those  frequencies  outside  lower  and  upper  cutoff  frequencies  (IAHR 
1987)  because  of  low  signal-to-noise  ratios.  Also,  RAO's  greater  than  10.0  or 
less  than  0.1  are  set  to  these  respective  upper  and  lower  limits. 

46.  The  stroke  control  signals  sc  for  each  of  the  61  paddles  which 
produced  the  measured  spectra  are  fourier  transformed  to  the  frequency  domain 
for  correction.  The  real  and  imaginary  fourier  coefficients  U  and  V  are 
divided  at  each  frequency  mAf  by  the  appropriate  value  of  the  smoothed  and 
averaged  RAO.  The  RAO  for  an  individual  gage  or  an  average  of  selected  gages 
can  be  used.  The  fourier  transform  relationship  is 
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s  (nAt) 
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exp 


/  -i2Trmn\ 

'  N  ' 


U(m)  -  iV(m)  =  At 


(33) 


* 


The  corrected  fourier  coefficients  are  then  inverse  fourier  transformed  to  the 
time  domain  to  give  the  new  control  signal 

N-1 

sc(nAt)  =  Af  [ U(m)  -  i V ( m )  ]  exp  (  (34) 

m=0 

The  square  root  in  Equation  32  is  necessary  because  spectral  estimates  are 
squared  quantities  of  the  fourier  coefficients 

S..(m)  =  ^  [u2(m)  +  V2(m)]  i  =  1,2,3, ...N  (35) 

r 

The  maximum  crest  and  minimum  trough  digital  values  are  2,048  and  -2,047, 
respectively.  The  digital  control  signals  for  the  DSWG  are  low-pass  filtered 
and  converted  to  analog  at  run  time. 

Wave  Grouping  Analysis 

47.  High  sea  waves  tend  to  fall  into  groups  rather  than  appear  individ¬ 
ually.  The  coastal  engineering  research  community  has  recognized  the  need  to 
model  wave  groups  as  well  as  spectral  waves.  High  waves  in  groups  can  produce 
more  damage  than  isolated  high  waves.  Engineers  are  finding  that  this  group¬ 
ing  is  important  in  the  motions  and  resonances  of  moored  structures  and  ves¬ 
sels,  harbor  resonance,  stability  and  overtopping  of  shore  protection  struc¬ 
tures,  and  surf  beat.  Because  of  the  nature  of  wave  grouping,  its  prediction, 
control,  and  analysis  are  especially  important  in  shallow-water  laboratory 
basins  such  as  CERC's  directional  spectral  wave  basin. 

48.  A  wave  group  is  defined  by  the  number  of  waves  which  exceed  some 
arbitrary  threshold,  such  as  the  mean  or  significant  wave  height.  Figure  2 
illustrates  the  concepts  of  a  run  of  high  waves  and  a  trtal  run.  The  length 
of  the  run  of  high  waves  is  the  number  of  waves  above  this  threshold.  The 
total  run  is  the  combination  of  the  run  of  high  waves  followed  by  the  run  of 
low  waves  (i.e.  succession  of  waves  which  fall  below  the  threshold  value). 

The  total  run  is  analogous  to  the  zero-crossing  period  of  the  wave  profile, 
except  that  the  time  series  is  composed  of  individual  wave  heights  rather  than 
surface  elevations.  Unless  specified  otherwise,  it  is  assumed  that  a  run  of 
high  waves  is  intended  when  a  run  length  s  mentioned. 


1  u  _ 


c- _ 


.•  _ 


TOTAL  RUNLENGTH 


ORDER  OF  WAVES 

Figure  2.  Wave  grouping  concepts 

predicting  wave  group  statistics  from  measured  wave  data  in  a  laboratory 
basin.  The  first  model  is  based  on  the  work  of  Goda  (1970).  It  assumes  that 
successive  wave  heights  are  uncorrelated  or  independent,  an  appropriate 
assumption  for  broadband  spectra.  Goda's  model  is  based  on  the  Rayleigh 
distribution  of  wave  heights,  which  assumes  a  narrow  band  spectrum.  Thus, 
this  model  suffers  from  an  internal  inconsistency  in  that  a  wave  field  cannot 
be  both  broad  and  narrow  banded.  While  Goda  (1976)  found  good  agreement  for 
broad  spectra  and  Elgar,  Guza,  and  Seymour  (19810  for  spectra  with  Qp  <  2  , 
Goda's  model  tends  to  break  down  for  narrow  spectra  as  the  run  length  remains 
constant . 

50.  Kimura's  model,  on  the  other  hand,  assumes  that  successive  wave 
heights  are  mutually  correlated  and  form  a  Markov  chain  (Kimura  1980).  The 
concept  of  mutual  correlation  implies  that  successive  waves  are  dependent  or 
correlated.  Kimura's  assumption  is  closer  to  what  is  observed  in  nature  and 
has  been  found  to  better  predict  observed  wave  grouping  due  to  sea  and  swell 
conditions  than  Goda's  model  (van  Vledder  1983a, b;  Thomas,  Baba,  and  Harish 
1986).  Successive  wave  heights  are  indeed  a  dependent  phenomenon.  A  high 
wave  rarely  appears  by  itself;  it  is  more  likely  to  be  followed  by  other  high 
waves.  It  is  as  though  the  waves  have  a  "memory,"  which  dictates  that  one 
high  wave  will  be  followed  by  another  high  wave  rather  than  a  low  wave. 

51.  Tn  the  paragraphs  below,  descriptions  of  the  probability  distribu¬ 
tions  and  run  lengths  for  a  run  of  high  waves  and  a  total  run  for  both  the 


Goda  and  Kimura  models  are  given.  The  interested  reader  is  referred  to  a 
report  and  paper  by  Briggs  ( 1987a, b)  for  more  complete  details. 

Goda 's  model 

52.  If  the  probability  p  that  a  wave  height  H  will  be  greater  than 
a  threshold  value  Hc  is 

p  =  Prob  [H  >  Hc]  (36) 

then  the  probability  that  this  will  not  occur  q  is  given  by 


q  =  Prob  [H  <  Hc]  =  1  -  p  (37) 

since  p  +  q  =  1  .  The  probability  distribution  (Goda  1985)  of  a  run  of 
successive  high  waves  P(J-j)  follows  the  geometric  distribution 

P(Jl)  =  pJl_1q  J,  =  1,2,3,...  (38) 


meaning  that  j^-1  successive  waves  will  exceed  the  threshold  and  the  J^th 
wave  will  not.  The  probability  of  successive  wave  heights  exceeding  the 
threshold  is  equal  to  the  product  of  separate  probabilities  for  each  event 
because  the  wave  heights  are  independent.  The  mean  run  length  of  a  run  of 
high  waves  J1  i3  given  by 


(39) 


tion 


and 


53.  For  a  total  run  of  J2  successive  waves,  the  probability  distribu- 
P(J2)  and  mean  run  length  j  are,  respectively 


Kimura1 s  model 

5^.  A  fundamental  assumption  of  Kimura' s  model  is  that  successive  wave 
heights  form  a  Markov  chain.  If  a  threshold  wave  height  Hc  is  selected, 
waves  with  height  H  will  fall  into  one  of  two  states:  state  1  where 


H  <  or  state  2  where  H  >  Hc  .  A  run  of  high  waves  begins  when  state  2  is 
first  reached.  A  transition  probability  matrix  [p]  is  given  by 


P11  p12 

[p]  =  (M2) 

-P21  P22. 


where  the  individual  elements  are  conditional  probabilities.  The  first  ele¬ 
ment  Pii  is  the  probability  that  neither  successive  wave  exceeds  the  thresh¬ 
old  H0  ,  and  p 22  is  the  probability  of  simultaneous  exceedance  by  both  wave 
heights. 

55.  These  conditional  probabilities  p^  and  1^22  are  defined  i-n 
terms  of  the  Rayleigh  distribution.  The  Rayleigh  probability  density  function 
(PDF)  for  individual  wave  heights,  and  the  joint  Rayleigh  PDF  for  two  succes¬ 
sive  wave  heights  are  calculated  numerically.  The  latter  requires  input  of  a 
correlation  parameter  <  which  is  related  to  the  correlation  coefficient 

Rhh  ,  a  measure  or  the  degree  of  correlation  or  dependence  between  wave 
heights. 

56.  The  correlation  coefficient  is  calculated  from  the  zero-meaned  mea¬ 
sured  wave  height  time  series.  The  autocorrelation  function  estimate  is  cal¬ 
culated  and  normalized  by  the  variance  to  give  the  correlation  coefficient 


Rhh(k)  =  2  N  -  k 

o 


H.  H.  , 
1  l+k 


k  =  1,2,3, 


where 

N  =  total  number  of  points  in  the  wave  height  time  series 
o  =  standard  deviation  of  wave  height  time  series 


The  lag  k  is  the  difference  in  number  between  wave  heights  and  is  equal 
to  one  for  successive  wave  heights.  For  every  other  wave  height,  the  corre¬ 
lation  coefficient  would  be  written  as  Rhh(2)  ,  every  third  wave  height 
Rhh(3)  ,  etc.  The  dependency  between  wave  heights  has  been  found  by  several 
investigators  (van  Vledder  1983a,  Arhan  and  Ezraty  1978,  Goda  1983)  to  fall 
off  rapidly  a3  the  lag  is  increased  beyond  successive  wave  heights  ( i . e . , 

57.  The  range  of  R^^( 1 )  is  0  to  1.0.  A  value  of  zero  corresponds  to 


the  Goda  model.  Several  investigators  (Goda  1985;  Thomas,  Baba,  and  Harish 
1986)  have  calculated  correlation  coefficients  of  0.24  to  0.3  for  successive 
wind  waves  and  0.5  to  0.8  for  swell.  The  amount  of  correlation  tends  to 
increase  with  higher  wave  heights  and  narrower  wave  spectra. 

58.  Once  R^(1)  *s  obtai-ne8  from  the  measured  data,  ie  is  calculated 
using  a  series  approximation  proposed  by  Battjes  (1974)  and  used  in  the  tran¬ 
sition  probabilities  p^  and  P22  •  This  approximation  is  very  good  for 
correlation  coefficients  less  than  0.7  to  0.8.  From  0.8  to  1.0,  the  differ¬ 
ence  between  the  approximation  and  the  exact  value  is  slight. 

59.  Thus,  for  Kiraura's  model  the  probability  of  a  run  of  successive 
high  waves  of  run  length  j1  is  defined  in  terms  of  the  transition  prob¬ 
ability  p22  as 

Jr1 

P(Jt)  =  P22  (1  -  P22)  J,  =  1,2,3 .  (44) 


The  mean  J1  is  similarly  defined  as 


r .  — ^ 

ji  1  -  p 


22 


(45) 


Note  the  similarity  of  Equations  44  and  45  with  Equations  38  and  39  for  Goda's 
model:  P22  and  (1  -  p22)  replace  p  and  q  ,  respectively. 

60.  For  a  total  run  of  Jg  successive  waves,  the  probability  distribu¬ 
tion  P(j2)  and  mean  run  length  are,  respectively, 


P(J2>  = 


(1  -  Pn)  (1  -  P22) 
P11  '  p22 


V1 

J11 


-  p 


V1 

22 


J2  =  1,2,3, 


(46) 


>2  =  '  -  P11  *  ’  -  F22 


(47) 


Again,  note  the  similarity  with  the  total  run  statistics  defined  for  Goda's 
model  in  Equations  40  and  41. 


PART  IV:  TEST  SETUP 


Wave  Basin  and  Generator 

61.  Spectral  waves  were  generated  in  CERC's  96-ft-long  by  121-ft-wide 
directional  spectral  wave  basin  (Figure  3).  It  can  accommodate  water  depths 
up  to  2  ft.  The  bathymetry  is  uniform  with  maximum  variations  of  0.35  and 
1.08  in.  over  a  5-ft  distance  and  the  entire  basin,  respectively  (Briggs  and 
Hampton  1987).  The  basin  has  a  unique  beach/wave  absorption  system,  the  pe¬ 
rimeter  of  which  is  lined  with  wave  absorber  frames  consisting  of  2-in.  layers 
of  rubberized  horsehair  installed  between  two  layers  of  expanded  sheet  metal. 
Each  segment  of  the  wave  absorber  beach  opposite  the  wave  generator  is  angled 
to  deflect  any  reflected  energy  away  from  the  primary  measurement  areas  toward 
the  sidewall  absorbers  (Figure  4).  Flume  tests  by  Briggs  and  Barnes  (sub¬ 
mitted)  for  a  range  of  wave  periods  and  heights  demonstrated  the  effectiveness 
of  the  frame  design.  The  average  measured  reflection  coefficient  of  12  per¬ 
cent  is  equivalent  to  a  comparable  gravel  beach  with  a  slope  of  6  to  1 1  deg 


Figure  3.  Schematic  of  basin  and  gage  arrangement 


Figure  4.  Photograph  of  wave  absorber  frames  along  beach 

which  could  be  installed  in  the  basin.  Additional  wave  absorption  is  provided 
behind  the  beach  segments  by  two  rows  (4  ft)  of  horsehair  extending  through 
the  water  column. 

62.  The  DSWG  was  designed  and  constructed  by  MTS  Systems  Corporation, 
Minneapolis,  Minnesota.  It  is  90  ft  long  and  consists  of  60  paddles  in  four 
modules  of  15  paddles  per  module,  each  1.5  ft  wide  and  2.5  ft  high.  It  is  an 
electronically  controlled  electromechanical  system.  The  paddles  are  indepen¬ 
dently  driven  at  each  of  the  61  joints  in  translational  motion  by  a  0.75-hp 
electric  motor.  Each  is  connected  with  a  flexible-plate  seal  to  provide 
continuity  and  minimize  the  introduction  of  spurious  waves  (Outlaw  and  Briggs 
1986). 

63.  The  modular  design  allows  one  or  more  of  the  modules  to  be  relo¬ 
cated  to  other  physical  model  test  sites  to  optimize  usage.  Because  of  this 
portability,  it  has  no  bottom  or  end  seals  and  is  referred  to  as  a  "wetback" 
design.  During  this  test  program  only  three  of  tne  modules  were  used. 

Thus,  the  total  length  of  the  DSWG  was  67.5  ft,  and  the  center  line  was  at 

Y  =  33.75  ft. 


Wave  Elevation  Measurement  System 


Gage  description 

64.  Seven  wave  gages  were  used  to  measure  surface  wave  elevations  in 
the  basin.  Initially,  parallel-wire  resistance  gages  were  used.  Because 
these  gages  are  sensitive  to  electronic  and  temperature  drift  and  cleanliness 
of  the  water,  they  were  replaced  with  capacitance  gages.  Flume  and  basin 
tests  were  conducted  prior  to  this  replacement  to  verify  the  accuracy  of  the 
capacitance  gages. 

Gage  arrangement 

65.  The  gages  were  mounted  in  aluminum  frames  to  minimize  inteference 
due  to  support  legs  and  arranged  in  a  directional  array.  Embedded  in  this 
array  was  a  2-3- 1 -7  linear  array  (Figure  3)  patterned  after  the  larger  linear 
array  design  at  CERC's  Field  Research  Facility  (FRF).*  The  linear  array  pro¬ 
vides  superior  resolution  capability  for  wave  components  at  or  near  the  same 
frequency  and  slightly  different  directions.  This  resolution  can  come  about 
if  the  waves  are  directional  or  have  been  refracted. 

66.  A  well-designed  linear  array  must  have  a  total  length  equal  to  the 
largest  wavelength  (i.e.  lowest  frequency)  at  the  largest  angle  expected.  It 
also  must  be  short  enough  to  avoid  aliasing  the  higher  frequency,  smaller 
wavelength  components.  Spatial  aliasing  occurs  when  half  the  wavelength  of 
the  highest  frequency  wave  does  not  exceed  the  distance  between  sensors.  It 
is  impossible  to  accurately  discern  the  smaller  wave  from  the  longer  wave. 
Thus,  there  is  a  classic  tradeoff  in  that  we  want  the  longest  array  to  opti¬ 
mize  resolution  of  the  low  frequency  components  and  the  shortest  spacing  to 
minimize  aliasing  of  the  high  frequency  waves.  The  secret  to  this  tradeoff  is 
to  select  a  minimum  distance  between  two  sensors  which  minimizes  aliasing  at 
the  high  frequency  cutoff  desired  while  simultaneously  providing  an  overall 
length  which  optimizes  resolution  at  the  low  frequency  cutoff.  By  clever 
arrangement  of  a  limited  number  of  gages,  both  objectives  can  be  achieved. 
Rather  than  spacing  the  gages  a  uniform  distance  apart,  they  are  spaced  at 
multiples  of  the  unit  lag  length  selected  above.  Thus,  all  wavelengths  be¬ 
tween  the  smallest  and  longest  are  covered  by  combinations  of  different  gages. 


*  From  hardnijt  at  seminar,  titled  "Linear  Arrays:  Wind-Wave  Directional 
Measurement  Systems,"  presented  by  J.  Oltman-Shay,  March  1987,  at  the  US 


For  the  2-3- 1 -7  array,  wave  periods  with  half  wavelengths  equal  to  1,  2,  3,  4, 
5,  6,  7,  8,  11,  and  13  lags  are  discernible. 

67.  The  procedure  consists  of  calculating  the  depth-limited  wavelength 
for  the  low  and  high  frequency  cutoffs  desired.  The  y-axis  wave  number  compo¬ 
nent  ky  is  calculated  as 


k  =  k  sin  9 

y 


(48) 


where  9  >  0  and  a  maximum  wave  direction  to  the  linear  array  of  ©  =  30  deg 
was  assumed.  The  corresponding  y-axis  wavelength  component  Ly  is  then 


(49) 


Finally,  the  array  length  capable  of  resolving  this  wavelength  is  equal  to 
one-third  to  one-half  of  Ly  .  Based  on  this  procedure,  a  lag  spacing  of 
16  in.  (1.33  ft)  was  selected  for  this  study. 

68.  The  two  remaining  gages  were  positioned  in  front  of  and  behind  the 
linear  array  for  redundancy  and  possible  reflection  analysis.  Table  1  sum¬ 
marizes  the  x-  and  y-axis  coordinates  of  each  of  the  seven  gages  in  the  direc¬ 
tional  array. 


Table  1 


Directional  Array  Coordinates 


Gage  Mo. 

X-axis,  ft 

Y-axis, 

1 

18 

33.75 

2 

20 

43.08 

3 

20 

33-75 

4 

20 

32.42 

5 

20 

28.42 

6 

20 

25.75 

7 

23 

33-75 

Gage  and  depth  calibration 


69.  The  measurement  rods  on  the  gages  had  an  overall  length  <->f  10  in. 


They  were  calibrated  over  a  6-in.  range  each  day  prior  to  conducting  tests. 6 


Jordan  controller  stepper  motor  was  used  to  automatically  raise  and  lower  the 
rod  through  a  series  of  11  steps  to  obtain  calibration  coefficients  using 
a  least  squares  linear  or  quadratic  fit.  This  averaging  technique,  using 
21  voltage  samples  per  gage,  minimi2es  the  effects  of  slack  in  the  gear  drives 
and  hysteresis  in  the  sensors.  Table  2  lists  the  quadratic  fit  and  maximum 
deviation  calibration  coefficients  (in  units  of  feet  times  10"-’)  for  each  gage 
for  each  day  of  testing.  The  water  depth  was  maintained  within  ±0.001  ft  of 
the  desired  1-ft  level  by  an  automatic  water  level  float  and  solenoid  control 
valve. 


Table  2 

Quadratic  Fit  Calibration  Coefficients  for  Capacitance  Wave  Gages 
Unidirectional  Spectral  Performance  Tests 

_ Maximum  Deviation,  10"^  ft _ 


Gage  No. 


No. 

Date 

RD0 1 

RD02 

RD03 

RD04 

RD05 

RD06 

RD07 

Min 

Max 

Ave 

1 

3/5/87 

314 

355 

79 

123 

390 

150 

367 

79 

390 

254 

2 

3/11/87 

81 

249 

75 

143 

97 

123 

274 

75 

274 

149 

3 

4/9/87 

Mnh 

4 

4/17/87 

5 

4/27/87 

121 

465 

357 

135 

198 

164 

417 

121 

465 

265 

Averages 

172 

356 

170 

134 

228 

146 

353 

92 

376 

223 

Data  Acquisition  and  Control  System 

70.  An  automated  data  acquisition  and  control  system  (ADACS)  is  used  to 
create  wave  generator  control  signals  for  each  of  the  60  wave  paddles,  monitor 
wave  paddle  displacement  feedback,  and  collect  and  analyze  time  series  experi¬ 
mental  data  from  the  wave  basin.  The  ADACS  is  controlled  by  a  DEC  VAX  11/750 
central  processing  unit  with  capability  i.  digital  control  signal  output 
through  an  IEEE  488  interface  and  120  channels  of  analog  to  digital  input. 

The  DSWG  control  signal  is  updated  at  a  rate  of  20  Hz.  Additional  information 
on  the  wave  basin,  wave  generating  system,  and  wave  gages  can  be  found  in  the 
report  by  Briggs  and  Hampton  (1987). 


Test  Case  Description 


71.  Thirty-four  test  cases,  characterizing  a  range  of  unidirectional 
spectral  parameters,  were  generated  to  quantify  the  performance  characteris¬ 
tics  of  the  DSWG.  These  cases  consisted  of  combinations  of  three  wave  peri¬ 
ods,  two  zero-moment  wave  heights,  three  spectral  peakedness  parameters,  and 
two  peak  wave  directions.  Wave  periods  were  1,  1.5,  and  2  sec,  wave  heights  1 
and  3  in.,  peakedness  parameters  3.3  and  7,  and  wave  directions  0  and  -30  deg. 
In  addition,  two  tests  with  peakedness  parameters  of  15  were  run  at  a  peak 
wave  direction  of  0  deg  to  investigate  the  effect  of  spectral  peakedness. 

Four  tests  with  wave  periods  of  1.0  sec  were  repeated  to  verify  repeatability. 
Finally,  four  tests  at  a  wave  period  of  1.0  sec  were  run  at  a  peak  direction 
of  30  deg  to  quantify  any  differences  in  basin  response  due  to  wave  orienta¬ 
tion.  Table  3  summarizes  the  spectral  parameters  for  these  34  test  cases. 

Test  Case  Parameters 


Simulation  and  generation  phase 

72.  The  D/A  rate  for  the  DSWG  is  20  Hz  (i.e.  time  increment  of 

0.05  sec),  so  61  time  series  of  10,000  points  or  500  sec  were  generated.  This 
control  signal  duration  corresponds  to  an  even  frequency  increment  of  0.002  Hz 
or  951  frequencies  between  the  lower  and  upper  cutoff  frequencies  of  0.1  t^ 

2.0  Hz.  The  DSA  amplitude  procedure  was  used  for  control  signal  generation. 
Data  collection  and  analysis  phases 

73.  After  a  waiting  time  of  25  sec  to  allow  slower  traveling  high  fre¬ 
quency  components  to  travel  beyond  the  farthest  gage,  wave  elevation  data  were 
sampled  at  10  Hz  (i.e.,  time  increment  of  0.1  sec).  A  minimum  of  200  waves  at 
the  peak  period  was  collected  as  recommended  by  Goda  (1985). 

74.  As  recommended  by  IAHR  (1986),  zero  downcrossing  analysis  was  used 
to  calculate  time  domain  parameters.  For  spectral  analysis,  the  data  were 
zero  meaned,  tapered  by  a  10  percent  cosine  bell  window,  and  band  averaged 
within  the  lower  and  upper  cutoff  frequencies  of  0.1  and  2.0  Hz,  respectively. 
The  resolution  bandwidth  is  a  function  of  peak  wave  period.  Every  10  spectral 
lines  were  averaged  to  give  20  degrees  of  freedom.  For  the  1.0-sec  peak  pe¬ 
riod,  the  resolution  bandwidth  was  0.05  Hz.  Table  4  lists  the  frequency 


Table  3 

Summary  of  Test  Cases 


Test 

Case 

Period 

Height 

Angle 

Gage 

ID 

sec 

in. 

Alpha 

Gamma 

deg 

Type* 

0-deg  Series 

1 

1 .0 

1 

0.00332 

3.3 

0 

R 

2 

1 .0 

1 

0.00227 

7.0 

0 

R 

3 

1.0 

3 

0.02990 

3.3 

0 

R 

4 

1.0 

3 

0.02070 

7.0 

0 

R 

5 

1.5 

1 

0.00113 

3.3 

0 

R 

6 

1.5 

1 

0.00084 

7.0 

0 

R 

7 

1.5 

3 

0.01030 

3.3 

0 

R 

8 

1.5 

3 

0.00770 

7.0 

0 

R 

9 

2.0 

1 

0.00059 

3.3 

0 

R 

10 

2.0 

1 

0.00045 

7.0 

0 

R 

1 1 

2.0 

3 

0.00538 

3.3 

0 

R 

12 

2.0 

3 

0.00410 

7.0 

0 

R 

-30  deg  Series 

13 

1.0 

1 

0.00332 

3.3 

-30 

C 

1  )l 

i  *■« 

1 .0 

1 

0.00227 

7.0 

-30 

C 

15 

1 .0 

3 

0.02990 

3.3 

-30 

C 

16 

i.O 

3 

0.02070 

7.0 

-30 

C 

17 

1  .5 

1 

0.00113 

3.3 

-30 

C 

18 

1.5 

1 

0.00084 

7.0 

-30 

C 

19 

1.5 

3 

0.01030 

3.3 

-30 

C 

20 

1.5 

3 

0.00770 

7.0 

-30 

C 

21 

2.0 

1 

0.00059 

3.3 

-30 

C 

22 

2.0 

1 

0.00045 

7.0 

C 

23 

2.0 

3 

0.00538 

3.3 

-30 

C 

24 

2.0 

3 

0.00410 

7.0 

-30 

C 

Gamma  15  Series 

25 

1.0 

1 

0.00144 

15.0 

0 

C 

26 

1 .0 

3 

0.01310 

15.0 

0 

c 

Repeat  Series 

27 

1.0 

1 

0.00332 

3.3 

0 

c 

28 

1.0 

1 

0.00227 

7.0 

0 

c 

29 

1 .0 

3 

0.02990 

3.3 

0 

30 

1.0 

3 

0.02070 

7.0 

0 

c 

Orientation  Series 

31 

1 .0 

1 

0.00332 

3.3 

30 

R 

32 

1 .0 

1 

0.00227 

7.0 

30 

R 

33 

1 .0 

3 

0.02990 

3.3 

30 

R 

34 

1.0 

3 

0.02070 

7.0 

30 

R 

«  R  = 

resistance 

gage,  C 

^  capacitance 

gage. 

Table  4 

Collection  and  Analysis  Parameters 


Description 


Data  Collection 


Peak  Period,  sec 


LP  =  1.5 


Record  length,  sec 
No.  points  in  record 
Frequency  increment,  Hz 


200 

2,000 

0.0050 


300 

3,000 

0.0033 


400 

4,000 

0.0025 


Spectral  Analysis 

Smoothed  bandwidth,  Hz  0.0500  0.0333 

No.  frequency  components  analyzed  481  721 

Unidirectional  Spectral  Analysis 


0.0250 

961 


No.  frequency  components  smoothed  23  35  47 


dependent  data  collection  and  spectral  analysis  parameters  for  each  of  the 
three  peak  wave  periods. 

75.  For  the  unidirectional  spectral  analysis,  the  data  were  again  zero 
meaned,  tapered  by  a  10  percent  cosine  bell  window,  and  smoothed  with  a 
Gaussian  smoothing  function  with  an  effective  width  of  0.05  Hz  between  the 
same  lower  and  upper  cutoff  frequencies.  Table  3  lists  the  number  of  fre¬ 
quencies  Gaussian  smoothed.  All  other  parameters  are  the  same  as  for  spectral 
analysis.  For  the  spreading  function  estimates,  a  resolution  of  5  deg  (i.e. 

72  increments  in  360  deg)  was  used.  The  predicted  values  were  calculated 
without  tapering  the  control  signals. 


PART  V:  TEST  RESULTS  AND  ANALYSIS 


Time  Series  Results 


Time  series  plots 

76.  The  measured  wave  elevation  time  series  plot  for  gage  7,  case  4  is 
shown  in  Figure  5.  Gage  7  is  representative  of  all  gages.  Plots  for  the 


other  33  cases  are  included  in  Appendix  B  (Volume  II).  The  x-axis  length  is 
200,  300,  or  400  sec  for  peak  period  cases  of  1.0,  1.5,  and  2.0  sec,  respec 


tively.  The  y-axis  scales  have  increments  of  0.05  and  0.10  ft  for  wave  height 
cases  of  1  and  3  in.,  respectively.  Comparisons  among  the  34  cases  illustrate 


the  effect  of  the  period,  height,  spectral  peakedness,  and  directional 


parameters . 


0.0  50.  100.  150.  200 


Time  (sec) 

CASE  4,  Tp  =  1 .0,  Hm0=3,  G  =  7.0,  D=0 
Figure  5.  Measured  wave  elevation  time  series 


Zero  downcros sing 

77.  Results  from  the  zero  downcrossing  analysis  for  each  case  are 
summarized  in  Table  5  which  is  organized  like  Table  3.  The  target  wave  period 
and  height  are  listed  for  reference.  The  measured  average  values  for  all 
seven  gages  of  the  significant  wave  period  T^/j  d  and  significant  wave 


Table  5 

Zero  Downcrossing  Analysis  Results 


Test 

Case 

ID 

Target 

Measured 

Average 

Normalized  Average 
Period  Height 

Period 

sec 

Height 

in. 

Period 

sec 

Height 
in . 

0-deg  Series 

1 

1.0 

1 .0 

0.88 

1 .07 

0.88 

1 .07 

2 

1.0 

1 .0 

0.91 

1 .04 

0.91 

1.04 

3 

1.0 

3.0 

0.91 

2.85 

0.91 

0.95 

4 

1 .0 

3.0 

0.93 

2.88 

0.93 

0.96 

5 

1.5 

1.0 

1.06 

0.83 

0.70 

0.83 

6 

1.5 

1.0 

1 . 16 

0.80 

0  77 

0.80 

7 

1.5 

3.0 

1.11 

2.61 

0.74 

0.87 

8 

1.5 

3.0 

1.18 

2.53 

0.78 

0.84 

9 

2.0 

1 .0 

1.29 

0.73 

0.64 

0.73 

10 

2.0 

1.0 

1.44 

0.70 

0.72 

0.70 

1 1 

2.0 

3.0 

1.31 

2.37 

0.65 

0.78 

12 

2.0 

3.C 

1 .44 

2.30 

0.72 

0.78 

-30  deg  Series 

13 

1.0 

1.0 

0.93 

0.82 

0.93 

0.82 

14 

1 .0 

1.0 

0.94 

0.84 

0.94 

0.84 

15 

1.0 

3.0 

0.93 

2.60 

0.93 

0.87 

’6 

1.0 

3.0 

0.96 

2.70 

0.96 

0.90 

17 

1  .5 

1 .0 

1 .21 

0.62 

0.81 

0.62 

18 

1  .5 

1 .0 

1.34 

0.67 

0.89 

0.67 

19 

1.5 

3.0 

1.28 

2.21 

0.85 

0.74 

20 

1.5 

3.0 

1.34 

2.23 

0.89 

0.75 

21 

2.0 

1 .0 

1.57 

0.61 

0.78 

0.61 

22 

2.0 

1.0 

1.70 

0.60 

n  Re 

v->  •  O  j 

0.60 

23 

2.0 

3.0 

1.59 

2.03 

0.^9 

0.68 

24 

2.0 

3.0 

1.72 

2.06 

0.86 

0.68 

Gamma  15  Series 

25 

1.0 

1 .0 

0.93 

0.98 

0.93 

0.98 

26 

1.0 

3.0 

0.94 

2.89 

0.94 

0.96 

Repeat  Series 

27 

1  .0 

1.0 

0.86 

1.03 

0.86 

1.03 

28 

1  .0 

1.0 

0.90 

1 .02 

0.90 

1.02 

29 

1.0 

3.0 

0.91 

2.83 

0.91 

0.94 

30 

1.0 

3.0 

0.94 

2.95 

0.94 

0.98 

Orientation 

Series 

31 

1.0 

1 .0 

0.90 

0.83 

O.^o 

0.83 

32 

1  .0 

1 .0 

0.93 

0.84 

0.93 

0.84 

33 

1  .0 

3.0 

0.^1 

2.43 

0.94 

0.81 

34 

1.0 

3.0 

0.94 

2.50 

0.91 

0.83 

two  columns.  It  is  the  ratio  of  the  measured  period  or  height  normalized  by 
the  target  value.  Appendix  C  lists  all  TSA  results  for  each  of  the  seven 
gages. 

78.  Measured  average  periods  are  plotted  against  target  period  for 
0-deg  (i.e.  1-12  and  25-26)  and  -30  deg  (i.e.  13-24)  cases  in  Figures  6a  and 
6b,  respectively.  The  cases  are  grouped  by  spectral  peakedness.  Normalized 
values  for  wave  period  ranged  from  a  low  of  0.64  for  case  9  at  2.0  sec  and 

0  deg  to  a  high  of  0.96  for  case  16  at  1.0  sec  and  -30  deg.  The  measured 
periods  for  the  eight  0-deg  cases  at  larger  periods  of  1.5  and  2.0  sec  seemed 
to  be  much  lower  than  their  target  values,  especially  relative  to  the  -30  deg 
series  cases.  The  normalized  period  values  for  these  eight  cases  were  0.11  to 
0.14  lower  than  their  counterparts  in  the  -30  deg  series.  Disregarding  these 
eight  cases,  the  lowest  normalized  period  was  0.78. 

79.  Measured  average  wave  heights  are  shown  in  Figures  7a  and  7b  for 
the  0-deg  and  -30  deg  series  cases,  respectively.  Again,  they  are  grouped  by 
gamma.  The  two  wave  height  levels  are  shown  individually.  As  previously 
observed  for  the  monochromatic  waves  (Briggs  and  Hampton  1987),  wave  height 
decreased  as  period  and  wave  direction  increased.  The  normalized  wave  heights 
ranged  from  a  low  of  0.60  for  case  22  at  2.0  sec  and  -30  deg  and  a  high  of 
1.07  for  case  1  at  1.0  sec  and  0  deg. 

80.  Results  from  the  repeatability  tests  are  shown  in  Figure  8.  It 
shows  the  measured  average  wave  heights  for  cases  1-4  (x-axis)  plotted  against 
the  corresponding  cases  27-30  (y-axis).  The  agreement  is  excellent.  Figure  9 
illustrates  the  effect  of  wave  orientation  on  the  measured  wave  height  for 
cases  13-16  at  -30  deg  (x-axis)  versus  cases  31-34  at  30  deg  (the  y-axis). 

The  agreement  is  very  good  with  a  slight  variation  for  the  3-in.  waves. 
Finally,  the  agreement  is  excellent  between  measured  wave  periods  for  both 
sets  of  cases. 

Spectral  analysis 

81.  Spectral  analysis  results  include  calculation  of  spectral  statis¬ 
tics  and  frequency  spectra  plots.  These  statistics  include  peak  wave  period 
Tp  ,  zero-moment  wave  height  H^q  ,  and  Coda  peakedness  parameter  Qp  . 

82.  Spectral  statistics.  Table  6  lists  statistics  for  the  spectral 

analysis  of  each  of  the  34  cases.  Again,  the  test  case  numbers  and  target 
wave  period  and  height  are  given  for  reference.  The  measured  average  of 
all  seven  gages  for  Tp  and  are  listed.  The  normalized  values  are 


MEASURED  PERIOD,  SEC  MEASURED  PERIOD,  SEC 


TARGET  PERIOD.  SEC 


11}  DOWNCROSSING  CASES:  CASES  1-12  4  25-26 
a.  Cases  of  0  deg 


0.5  1.0  1.5  2.0  2.5 

TARGET  PERIOD.  SEC 

12  DOWNCROSSING  CASES:  CASES  13-24 

b.  Cases  of  -30  deg 

Figure  6.  Measured  versus  target  zero-crossing  wave  periods 
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Figure  7.  Measured  versus  target  zero-crossing  wave  heights 
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Figure  8.  Repeatability  series  of 
zero-crossing  wave  heights 


MEASURED  HEIGHT.  IN 

4  DOWNCROSSING  CASE 5 :  CASES  13-16  VS.  31- 

Figure  9.  Orientation  series  of 
zero-crossing  wave  heights 


Table  6 

Spectral  Analysis  Results 


'est 

Target 

Measured  Average 

Goda 

lase 

Period 

Height 

Period 

Height 

Normal i2ed  Average 

Spectral 

in. 

sec 

in. 

Period 

Height 

Peakedness 

ft 

•deg  Series 

1 

1.0 

1.0 

1.15 

0.97 

1.15 

3. 18 

2 

1.0 

1.0 

0.99 

1.11 

0.99 

1.11 

4.11 

3 

1.0 

3.0 

1.01 

2.93 

1.01 

0.98 

3.24 

4 

1.0 

3.0 

2.96 

1 .02 

0.99 

4.23 

5 

1.5 

1.0 

1.44 

0.90 

0.96 

0.90 

1 .87 

6 

1.5 

1.0 

1.51 

0.87 

1.01 

0.87 

2.28 

7 

1.5 

3.0 

1.40 

2.76 

0.94 

0.92 

1.85 

8 

1.5 

3.0 

1.44 

2.68 

0.96 

0.89 

2.13 

9 

2.0 

1.0 

1.92 

0.80 

0.96 

0.80 

1.63 

10 

2.0 

1.0 

1.92 

0.77 

0.96 

0.77 

2.20 

11 

2.0 

3.0 

1.97 

2.54 

0.99 

0.85 

1.54 

12 

2.0 

3.0 

1.97 

2.48 

0.99 

0.83 

2.04 

-30  deg  Series 

13 

1.0 

1.0 

0.86 

1.00 

0.86 

4.37 

14 

1.0 

1.0 

0.99 

0.89 

0.99 

0.89 

5.83 

15 

1.0 

3.0 

0.98 

2.67 

0.98 

0.89 

4.09 

16 

1.0 

3.0 

1.02 

2.75 

1.02 

0.92 

5.38 

17 

1.5 

1.0 

1.46 

0.66 

0.97 

0.66 

2.34 

18 

1.5 

1.0 

1.51 

0.72 

1 .01 

0.72 

3.49 

19 

1.5 

3.0 

1.49 

2.32 

0.99 

0.77 

2.40 

20 

1.5 

3.0 

1.52 

2.35 

1.01 

0.78 

3. 14 

21 

2.0 

1.0 

1.94 

0.65 

0.97 

0.65 

2.21 

22 

2.0 

1.0 

1.94 

0.64 

0.97 

0.64 

3.13 

23 

2.0 

3.0 

1.96 

2.16 

0.98 

0.72 

2.17 

24 

2.0 

3.0 

1.98 

2.18 

0.99 

0.73 

3.19 

Gamma  15  Series 

25 

1.0 

1.0 

0.98 

1.06 

0.98 

1.06 

6.59 

26 

1.0 

3.0 

1.01 

2.96 

1.01 

0.99 

5.92 

Repeat  Series 

27 

1 .0 

1.0 

1.11 

1 .00 

1.11 

3.00 

28 

1.0 

1.0 

0.99 

1.09 

0.99 

1 .09 

3.88 

29 

1.0 

3.0 

0.99 

2.89 

0.99 

0.96 

3.32 

30 

1.0 

3.0 

1.02 

3-01 

1.02 

1.00 

4.41 

Orientation  Series 

31 

1 .0 

1.0 

1 .00 

0.88 

1.00 

0.88 

4.06 

32 

1.0 

1 .0 

0.99 

0.90 

0.99 

0.90 

6.11 

33 

1.0 

3.0 

1.00 

2.51 

1.00 

0.84 

3.76 

34 

1.0 

3.0 

1.01 

2.55 

1.01 

0.85 

5.45 

calculated  as  before.  Goda's  Qp  is  listed  in  the  last  column  for  use  in  the 
wave  grouping  analysis. 

83.  Measured  average  periods  are  plotted  against  target  periods  for 
0-dcg  (i.e.  1-12  and  25-^6;  and  -30  deg  (i.e.  13-24)  cases  in  Figures  10a  and 
10b,  respectively.  They  are  analogous  to  Figures  6a  and  6b.  The  range  of 
normalized  periods  is  from  a  low  of  0.94  for  case  7  at  1.5  sec  and  0  deg  to  a 
high  of  1.02  for  cases  4,  16,  and  30  at  1.0  sec.  This  variation  in  measured 
periods  is  much  less  than  the  corresponding  zero-crossing  values. 

84.  Measured  average  wave  heights  are  shown  for  the  0-deg  and  -30  deg 
series  cases  in  Figures  11a  and  11b,  respectively.  They  are  analogous  to 
Figures  7a  and  7b  for  the  zero-crossing  analysis.  The  cases  are  grouped  by 
gamma  and  wave  height.  The  normalized  heights  varied  from  a  low  of  0.64  for 
case  22  at  2.0  sec  and  -30  deg  to  a  high  of  1.15  for  case  1  at  1.0  sec  and 

0  deg.  Again,  wave  height  decreases  with  increasing  wave  period  and  direc¬ 
tion.  Even  though  these  values  differed  in  period  and  height  from  their 
target  values,  they  are  much  closer  than  the  corresponding  zero-crossing 
analysis  results.  The  RAO  transfer  functions  should  correct  the  control 
signals  to  give  the  required  wave  heights. 

85.  Figures  12  and  13  show  the  repeatability  and  wave  orientation  test 
series  wave  heights,  respectively.  The  format  is  the  same  as  for  Figures  8 
and  9  for  the  zero-crossing  analysis.  The  variations  in  wave  period  and 
height  for  the  repeat  and  wave  orientation  series  with  their  counterparts  were 
insignificant. 

86.  Frequency  spectra.  Figure  14  illustrates  the  measured  versus  pre¬ 
dicted  frequency  spectra  for  case  4:  Tp  =  1.0  sec,  =  3  in.,  y  =  7.0  , 

and  D  =  0  deg.  This  case  is  used  throughout  this  report  as  an  example.  Plots 
of  the  other  cases  are  included  in  Appendix  D.  The  x-axis  is  from  0  to 

2.0  Hz.  The  y-axis  varies  according  to  wave  height.  The  range  is  0  to  3  for 
1-in.  waves  and  0  to  24  for  3- in.  waves.  The  dashed  line  is  the  predicted 
spectrum  and  the  solid  line  gives  the  measured  spectrum.  The  agreement  is 
excellent. 

87.  The  agreement  for  the  1.0-sec  cases  was  generally  excellent  for  all 
cases  regardless  of  wave  height,  direction,  or  spectral  peakedness.  The  shape 
of  the  spectra  for  the  longer  period  cases  was  not  as  accurate  as  desired. 

The  measured  spectral  peak  in  these  cases  was  usually  less  than  the  predicted 
value  by  a  factor  of  the  order  of  2.  The  magnitude  of  the  low  frequency  tail 


moo 


4  SPECTRAL  CASES:  CASES  1-4  VS.  27-30 
Figure  12.  Repeatability  series  of  spectral  wave  heights 


4  SPECTRAL  CASES:  CASES  13-16  VS.  31-34 


Figure  13.  Orientation  series  of  spectral  wave  heights 


TP-1.00.  HM0=3. 00,  3-7.0,  D=0.  3W=.05 
Figure  14.  Measured  versus  predicted  frequency  spectra 

energy  is  relatively  small.  Any  deviations  from  the  predicted  values  can  be 
the  result  of  low  SNR's  ( signal-to-noise  ratio)  or  the  presence  of  bounded  or 
free  long  waves.  The  RAO's  (see  paragraph  93)  indicated  substantial  differ¬ 
ences  which  are  hidden  in  the  frequency  spectra  graphs  because  of  the  linear 
scale.  The  high  frequency  tail  energy  for  all  cases,  however,  was  never  more 
than  two  times  greater  than  the  predicted  value.  The  largest  variations 
occurred  for  an  of  1  in.  The  RAO  transfer  functions  should  be  able  to 

correct  for  the  observed  differences  in  one  or  more  iterations. 

Unidirectional  Spectral  Analysis 

88.  Unidirectional  spectral  analysis  results  include  peak  frequency 
estimates,  plots  of  directional  spreading  for  all  34  cases  at  the  peak  fre¬ 
quency,  and  a  comparison  of  the  DSA  and  NSA  methods  of  spectral  simulation  for 
two  cases. 

Peak  frequency  estimates 

89.  The  measured  peak  frequencies  agreed  with  the  target  values  in  all 
cases  within  +0.01  Hz. 


Directional  spreading  function 

90.  The  measured-versus-predicted  directional  spreading  function  esti¬ 
mates  for  case  f!  arc  ^I.own  in  figure  15.  The  spreading  functions  for  the 
other  33  cases  are  included  in  Appendix  E.  The  average  of  all  seven  gages  is 
used  for  the  measured  results  and  of  all  61  paddles  for  the  predicted  values. 
The  truncation  of  the  fourier  series  tends  to  produce  negative  side  lobes  as 
an  artifact  of  the  fourier  transform  method.  Thus,  negative  spiaading  func¬ 
tion  estimates  can  sometimes  be  produced  at  wave  directions  on  either  side  of 
the  mean  wave  direction.  These  negative  values  have  been  suppressed  and  set 
to  zero  on  these  plots. 
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CASE  4:  U 1 22 1 1 ,  FP=1.00,  HM0=3.  G=7.0.  0=0.  S=1 
Figure  15.  Measured  versus  predicted  directional  spreading 

91.  These  plots  indicate  the  amount  of  directional  spreading  and  the 
mean  wave  direction  at  the  peak  frequency.  In  all  cases  the  measured  value  of 
mean  wave  direction  agrees  with  the  predicted  values  of  0,  30,  or  -30  deg. 

The  agreement  of  the  measured  spreading  with  the  predicted  spreading  is  excel¬ 
lent  for  the  cases  with  1.0-sec  peak  period.  The  larger  period  cases  indicate 
a  wider  spreading  than  predicted,  which  increases  with  increasing  peak  period. 
A  possible  explanation  for  the  differences  could  be  a  frequency-dependent 
artifact  of  the  simulation  and  analysis  programs.  A  different  random  number 


seed  in  the  simulation  program  might  produce  a  different  realization  with  more 
stable  characteristics.  Also,  the  analysis  program  might  be  introducing  arti¬ 
ficial  smoothing  of  the  spreading  estimates.  MLM  and  Maximum  Entropy  Method 
(MEM)  are  analysis  techniques  which  tend  to  give  higher  resolution.  The  mea¬ 
sured  data  might  be  analyzed  using  both  of  these  methods  to  determine  if  more 
stable  results  can  be  obtained.  Another  possible  explanation  of  this  differ¬ 
ence  in  measured  spreading  may  be  the  departure  of  the  frequency  spectra  from 
the  desired  value.  Since  the  directional  wave  spectrum  is  the  product  of  the 
frequency  spectra  and  the  spreading  function,  any  change  in  the  frequency 
spectrum  will  cause  a  corresponding  change  in  the  spreading  function.  By 
increasing  the  spectral  peak  energy  using  the  appropriate  RAO,  the  spreading 
function  may  increase  and  become  narrower. 

PSA  and  NS A  comparison 

92.  A  comparison  of  the  DSA  and  NSA  amplitude  methods  for  simulating 
unidirectional  wave  spectra  was  made  for  two  cases.  Target  spectral  param¬ 
eters  were  Tp  =  1.0  sec  ,  =  2.0  in.,  y  =  7.0  ,  and  0m  =  0  deg.  This 

case  is  similar  to  case  4  except  that  the  wave  height  is  smaller.  It  corre¬ 
sponds  to  a  very  narrow  directional  spectrum  in  both  frequency  and  direction: 
conditions  most  necessitating  use  of  the  NSA  method.  Figures  16  and  17  show 
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CASES  35  4  36:  TP=1.0,  HM0=2.  G=7.0.  D=0.  S=1 
Figure  16.  Measured  frequency  spectra  for  DSA 
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CASES  35  &  36:  TP=1.  HM0=2,  6=7.0.  D=0,  S=1 

Figure  17.  Measured  spreading  functions  for  DSA 
and  NSA  methods 

the  measured  frequency  spectra  and  spreading  functions,  respectively,  for  both 
methods.  Even  though  the  frequency  spectrum  for  the  NSA  method  is  larger  than 
that  for  the  DSA,  the  measured  wave  height  for  each  was  2.04  and  2.16  in., 
respectively.  Both  values  are  statistically  acceptable.  The  NSA  value  repre¬ 
sents  a  larger  realization  within  the  realm  of  probability.  The  peak  wave 
periods  for  both  methods  agree  with  the  predicted  value.  Finally,  the  mea¬ 
sured  spreading  functions  show  little  variation  with  each  other.  Thus,  the 
DSA  method  is  a  satisfactory  method  for  simulating  directional  wave  spectra  in 
laboratory  basins. 


Response  Amplitude  Operator  Functions 


93.  Figure  18  is  a  semilog  plot  of  the  RAO  transfer  function  for 
case  4.  The  RAO's  for  the  other  cases  are  included  in  Appendix  F.  The  RAO's 
for  each  gage  and  the  average  for  all  seven  are  shown.  The  RAO  estimates  were 
Gaussian  smoothed  to  0.05  Hz,  and  all  estimates  within  contiguous  0.05-Hz 
bands  were  averaged.  Values  less  than  1.0  indicate  that  the  spectrum  has  to 
be  increased  at  a  nart.icnlar  frpnnpnrv  Accnrrlinul  u  uaUiA?  arestpr  t  ban  10 
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T  P  = 1 . 00 ,  HMQ  =  3 . 00 ,  G  =  7 . 0 ,  0  =  0.  BW= . 05 

Figure  18.  RAO  transfer  function 

indicate  the  opposite  trend.  As  discussed  previously,  the  RAO  will  be  used  in 
future  tests  to  adjust  the  generated  spectrum  to  compensate  for  DSWG  and  basin 
characteristics.  This  RAO  correction  will  be  especially  helpful  for  the  1.5- 
and  2.0-sec  peak  period  spectral  cases. 


Wave  Grouping  Analysis 


94 .  In  this  section,  comparisons  between  God3's  and  Kimura's  wave  group 
analysis  methods  for  a  run  of  high  waves  are  presented  and  discussed.  The 
wave  height  time  series  plots  from  the  zero-crossing  analysis  are  presented 
first.  Next,  plots  of  correlation  coefficients,  probability  distributions, 
mean  run  lengths,  total  run  length,  and  Goda's  peakedness  parameter  for  both 
methods  are  presented. 

Wave  height  time  series 

95.  Figure  19  is  a  plot  of  the  wave  height  time  series  for  case  4. 
Appendix  G  contains  plots  for  the  other  cases.  The  x-axis  is  individual  zero- 
crossings  or  waves  in  a  record  instead  of  time.  The  y-.,xis  is  scaled  with 

;.c  foments  of  n .  5  or  ’  .  f.)  in.  depending  on  the  target  wave  height  .  Rot  n  me  r 
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CHANNEL  NUMBER  =  7 

Figure  19-  Measured  wave  height  time  series 

H  (instead  of  in  these  plots)  and  zero-moment  wave  height  thres¬ 

holds  are  shown.  The  measured  number  of  crossings  above  these  thresholds  is 
used  in  the  wave  grouping  analysis. 

Correlation  coefficients 

96.  Figures  20a  to  20c  show  the  effect  of  various  lag  lengths  k 

between  successive  wave  heights  on  Rhh(k)  for  0-deg  cases  with  periods  of 
1.0,  1.5,  and  2.0  sec,  respectively.  Each  case  is  plotted  separately  in  each 
figure.  Figures  21a  to  21c  are  analogous  plots  for  the  -30  deg  cases.  The 
RhhOO  for  the  1.5-  and  2.0-sec  period  cases  at  0  deg  may  be  artificially  low 
because  of  the  low  directional  spectra  recorded.  As  expected,  the  correlation 
drops  off  rapidly  for  lags  greater  than  one  for  all  cases.  The  spectral  cases 
with  larger  y  values  have  a  higher  Rhh(k)  .  Calculated  values  of  Rhh(i) 

for  wind  seas  cases  (i.e.  y  =  3.3)  generally  agree  with  those  of  previous 

investigators  (see  paragraph  57).  The  cases  for  Tp  =  1.0  sec  are  high,  and 

those  for  Tp  =  1.5  and  2.0  sec  at  0  deg  are  low.  For  the  swell  cases 

(i.e.  y  =  7.0  and  15.0),  the  agreement  is  good  except  for  low  values  for  Tp 
=  1.5  and  2.0-sec  cases. 

97.  The  Rhh(k)  for  the  30-deg  cases  are  uniformly  higher  than  the 
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Figure  20.  Effect  of  lag  length  on  correlation  coefficient  at  0  deg 


corresponding  0-deg  values.  A  cursory  examination  of  the  measured  values  for 
Rhh^k)  indicates  a  trigonometric  relationship  between  0-  and  30-deg  cases. 
Also,  the  Tp  =  1.5-  and  2.0-sec  cases  show  an  upturn  in  R^(k)  for  lags 
greater  than  four  indicating  a  correlation  for  every  fifth  wave.  Further 
investigation  into  these  relationships  is  warranted. 

Probability  distributions 

98.  Figures  22a  and  22b  show  the  probability  for  a  run  of  high  waves 

exceeding  a  specified  run  length  for  thresholds  of  1^  and  >  respec¬ 

tively,  for  case  4.  The  predicted  Goda  and  Kimura  probabilities  are  also 
plotted  on  each  figure.  The  theoretical  Goda  line  is  a  constant  which  only 
depends  on  the  threshold  value  used.  The  theoretical  Kimura  line,  however,  is 
based  on  an  average  R^(k)  value  for  all  lags,  weighted  mainly  by  Rhtt(  1)  . 

A  value  of  Rhh  =  0.4  was  use:1  in  Figures  22a  and  22b.  The  calculated  prob¬ 
abilities  are  predicted  better  by  the  Kimura  model,  regardless  of  threshold 
value. 

99.  The  probability  distributions  for  the  other  cases  are  contained 
in  Appendix  H.  Values  of  Rhh  of  0.2,  0.3,  0.4,  0.5,  and  0.6  were  used  to 
calculate  the  theoretical  Kimura  probabilities  based  on  the  measured  Rhh(k) 
(see  paragraphs  96  and  97).  As  expected  (paragraph  49),  the  Goda  model  better 
predicted  the  broader  spectra  indicated  by  the  lower  Rhh(k)  values  of  0.2. 
Mean  run  length 

100.  Figures  23a  and  23b  show  the  relationship  between  the  mean  run 
length  Jj  for  a  threshold  of  1^  and  the  correlation  coefficient  R^(1) 
for  all  0-  and  30-deg  cases,  respectively.  Figures  24a  and  24b  show  the 
corresponding  relationship  for  a  threshold  of  .  The  Goda  and  Kimura 

theoretical  values  are  also  shown  in  each  figure.  The  measured  values  more 
closely  match  the  Kimura  predictions,  although  they  are  slightly  lower.  The 
difference  between  the  two  models  is  greater  for  the  threshold.  Goda's 

model  gives  a  constant  value  for  different  values  of  Rhh(D  that  seriously 
underpredicts  the  degree  of  wave  grouping  except  for  broad  spectral  cases 
with  R^(  <  0.2  .  These  results  are  in  agreement  with  those  found  by 

previous  investigators  (Goda  1983;  Rye  1974;  Kimura  1 980 ;  van  Vledder  1983a, b; 
Elgar,  Guza,  and  Seymour  1984;  Thomas,  Baba,  and  Harish  1986).  The  values 
of  Rhh(1)  and  J1  are  slightly  higher  for  the  30-deg  cases. 

Total  run  length 


101.  Figures  25a  and  25b  illustrate  the  relationship  between  total  run 
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Figure  22.  Probability  distribution  for  run  of  high  waves 


MEF.N  HUN  LENGTH,  H  >  HM  MERN  RUN  LENGTH, 


CCGGZLGT i GN  COEFFICIENT.  RHH  I  i 


a.  Cases  of  0  deg 


1  .  Li  ‘ - - - J— — - L— - - - ‘ - i- - * - < 

0.0  0.2  0.4  0.6  0.8  1  .  C 


OPRE^'-T  ’  ON  COEFR  i  C  I  ENT  .  dhh  ■.  \ 


b.  Cases  of  30  deg 


Figure  23.  Run  length  versus  correlation  coefficient, 
mean  height  threshold 
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length  J2  and  R^( 1)  for  the  0-  and  30-deg  cases,  respectively.  Values 
for  both  thresholds  are  shown  in  each  figure.  The  data  for  the  30-deg  cases 
substantiate  the  trends  of  the  0-deg  cases  For  a  threshold  of  ,  a  strong 
linear  or  quadratic  trend  is  indicated.  For  the  threshold,  the  data 

show  more  scatter.  A  least  squares  fit  could  easily  be  calculated  for  both 
thresholds . 

Goda's  peakedness  parameter 

102.  Goda  (1970,  1976)  proposed  his  spectral  peakedness  factor  Qp  as 
a  possible  means  of  quantifying  the  degree  of  wave  groupiness  present.  Its 
value  ranges  from  1  for  white  noise,  2  for  wind  waves,  and  4-8  for  swell  con¬ 
ditions.  Previous  investigators  (Yamaguchi  1981,  Kimura  1980),  in  addition  to 
Goda,  found  that  the  average  group  length  increases  as  Qp  increases,  and  a 
narrow  spectrum  has  a  greater  degree  of  grouping  than  a  wide-banded  spectrum. 
Long  (1987),  however,  found  that  Qp  could  have  the  same  value  for  many  dif¬ 
ferent  spectral  shapes.  Thus,  although  Qp  is  probably  not  the  best  indica¬ 
tor  of  wave  groupiness,  it  does  show  general  trends. 

103.  Figures  26a  and  26b  illustrate  the  relationship  between  Kimura' s 

predicted  and  Qp  fcr  the  0-  and  30-deg  cases,  respectively.  Both 

threshold  values  are  shown  in  each  figure.  The  general  trends  for  both  thres¬ 
holds  are  supported  for  both  0-  and  30-deg  cases.  A  strong  linear  or  quad¬ 
ratic  trend  is  evident  in  the  0-deg  cases  for  both  thresholds.  There  is  more 
scatter  for  Qp  >  4  for  the  30-deg  cases,  however.  The  data  substantiates 
the  trend  toward  an  increase  in  as  Qp  increases. 

104.  Finally,  Rhh( 1)  is  plotted  as  a  function  of  Qp  for  the  0-  and 
30-deg  cases  in  Figures  27a  and  27b,  respectively.  The  solid  line  in  each 
figure  is  the  observed  values  of  Goda  (1976).  The  measured  values  agree  with 
this  data  set.  An  increase  in  Qp  indicates  an  increase  in  Rhh( 1 )  ,  with  a 
possible  leveling  off  or  decrease  at  Qp  >  5.5  . 
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Figure  27.  Correlation  coefficient  versus  Goda's 
peakedness  parameter 


PART  VI:  SUMMARY  AND  RECOMMENDATIONS 


105.  The  DSWG  is  an  important  resource  for  the  study  of  naturally 
occurring  sea  conditions  in  laboratory  physical  models.  Coastal  structures 
respond  differently  to  unidirectional  spectral  waves  than  they  do  to  mono¬ 
chromatic  waves.  The  modular  design,  portability,  method  of  paddle  connec¬ 
tion  and  displacement,  and  electrical  power  of  the  DSWG  produce  clean,  quiet 
unidirectional  spectra  from  any  peak  direction  without  having  to  physically 
relocate  the  wave  maker.  Thus,  transformation  of  wave  energy  in  shallow  and 
intermediate  water  depths  due  to  shoaling,  refraction,  and  diffraction  can  be 
more  efficiently  studied. 

106.  A  series  of  36  unidirectional  spectral  waves  was  generated  to 
quantify  performance  characteristics  of  the  DSWG  basin.  This  report  deals 
with  the  mechanics  of  simulating,  generating,  measuring,  and  analyzing  unidi¬ 
rectional  wave  spectra.  This  assessment  is  a  necessary  first  step  in  our 
overall  research  program  to  ensure  accurate  generation  and  analysis  of  natu¬ 
rally  occurring  waves  prior  to  site-specific  tests  with  coastal  structures 
such  as  harbors,  breakwaters,  Jetties,  mounds,  etc.  The  test  series  consisted 
of  combinations  of  three  wave  periods,  two  zero-moment  wave  heights,  three 
spectral  peakedness  parameters  (narrowness  of  spectrum),  and  two  peak  wave 
directions  covering  a  range  of  wave  environments.  Wave  periods  were  1.0,  1.5, 
and  2.0  sec;  wave  heights  1  and  3  in.;  peakedness  parameters  3-3,  7.0,  and 
15.0;  and  peak  wave  directions  0  and  -30  deg.  In  addition,  four  tests  were 
repeated  to  verify  repeatability,  and  four  tests  were  run  at  a  peak  direction 
of  30  deg  to  investigate  the  effect  of  wave  orientation  in  the  basin.  Water 
depth  for  all  tests  was  1.0  ±  0.001  ft. 

107.  Measurements  were  made  at  seven  locations  within  the  basin  by 
resistance  and  capacitance  wave  gages.  Comparison  of  results  between  the  two 
types  of  gages  in  a  flume  and  in  the  basin  indicated  excellent  agreement.  The 
gages  were  arranged  in  a  2-3- 1 -7  linear  array  embedded  in  a  directional  array 
in  the  shape  of  a  cross.  It  was  centrally  located  20  ft  in  front  of  the  DSWG 
along  its  center  line.  A  unit  lag  spacing  of  1.33  ft  was  selected  to  optimize 
resolution  while  minimizing  spatial  aliasing. 

108.  Unidirectional  spectra  were  simulated  as  the  product  of  a  unimodal 
TMA  frequency  spectrum  and  an  empirically  derived  wrapped  normal  directional 
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frequency  domain  using  a  double  summation,  deterministic  amplitude,  random 
phase  model.  Another  method,  the  non-deterministic,  random  phase  model,  pro¬ 
duces  an  exact  Gaussian  realization,  but  at  the  cost  of  more  variability  and 
computer  time.  Two  cases  were  generated  to  compare  the  two  methods.  Results 
showed  little  difference  for  the  spectral  characteristics  simulated. 

109.  Unidirectional  spectral  analysis  was  based  on  the  fourier  series 
expansion  of  the  directional  spreading  function.  The  measured  surface  data 
for  each  gage  was  zero-meaned,  windowed,  fourier  transformed  with  a  "235"  FFT, 
and  Gaussian  smoothed  in  the  frequency  domain  to  obtain  the  cross-spectral 
matrix  of  autospectra  and  cross-spectra.  This  matrix  is  substituted  into  the 
parameterized  spreading  function  to  obtain  a  set  of  Ng  simultaneous,  linear 
equations.  A  least  squares  fourier  transform  method  for  numerical  integration 
is  used  to  invert  a  matrix  of  independent  variables  to  solve  for  the  fourier 
coefficients  of  the  spreading  function.  The  technique  for  doing  this  is  a 
linear,  stepwise  regression  model. 

110.  Because  of  leakage  around  and  under  a  "wetback"  wave  maker,  elec¬ 
tronic  and  mechanical  losses,  use  of  linear  wave  theory,  and  basin  response 
characteristics,  the  measured  spectra  may  not  faithfully  reproduce  the  target 
spectra.  The  height-to-stroke  transfer  function  only  compensates  for  an  ideal 
wave  maker.  Thus,  an  RAO  transfer  function  was  calculated  for  each  spectral 
control  signal  to  compensate  for  observed  variations  in  peak  period,  wave 
height,  or  spectral  shape.  The  RAO  is  the  ratio  of  measured-to-predicted 
spectral  density  for  each  gage  at  each  frequency.  A  Gaussian  smoothed,  aver¬ 
aged  RAO  value  will  be  used  to  correct  the  original  stroke  control  signal  in 
the  frequency  domain  prior  to  generating  the  spectra  again. 

111.  High  waves  in  groups  can  produce  more  damage  than  isolated  high 
waves.  Wave  grouping  is  important  in  the  motions  and  resonances  of  moored 
vessels,  harbor  resonance,  structure  stability,  and  overtopping  of  shore  pro¬ 
tection  structures.  Thus,  it  is  necessary  to  accurately  calculate  and  predict 
wave  group  statistics.  Both  Coda's  and  Kimura's  models  were  compared  using 
measured  wave  height  time  series.  Goda's  model  assumes  that  successive  wave 
heights  are  independent,  an  appropriate  assumption  for  broad-banded  spectra. 
Kimura's  model  assumes  mutual  correlation  or  dependence  of  successive  heights, 
a  condition  which  is  most  often  observed  in  nature. 

112.  A  minimum  of  200  waves  at  the  peak  period  was  collected  with  a 
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spectral,  frequency  response,  and  wave  grouping  analyses  were  performed  for 
each  case.  Averages  for  all  seven  gages  of  wave  period,  height,  direction, 
spectral  shape,  directional  spreading,  RAO  transfer  functions,  and  wave  group 
statistics  were  calculated. 

113.  Measured  values  of  zero-crossing  wave  period  th1/3  d  and 
H1/3  d  were  averaged  and  normalized  by  target  values.  The  zero-crossing 
values  were  lower  than  their  corresponding  opectral  analysis  values.  The 
normalized  values  of  TH1/^  ^  ranged  from  a  low  of  0.64  to  a  high  of  0.96. 

The  1.0-sec  wave  period  cases  agreed  well  with  the  target  values.  The  1.5- 
and  2.0-sec  period  cases  at  0  deg  were  the  lowest.  Normalized  values  of 
H1/3,d  rar|ged  from  a  low  of  0.60  to  a  high  of  1.07.  Wave  height  decreased 
with  increasing  wave  direction  and  period. 

114.  For  spectral  analysis,  the  data  were  zero-meaned,  windowed, 
fourier  transformed,  and  band  averaged  between  lower  and  upper  cutoff  frequen¬ 
cies  for  20  degrees  of  freedom.  The  measured  spectral  analysis  peak  period 

Tp  and  zero-moment  wave  height  Hm0  were  also  averaged  and  normalized  by 
target  values.  The  agreement  between  measured  and  target  Tp  was  excellent. 
Normalized  Tp  ranged  between  0.94  to  1.02.  The  normalized  ranged  from 

a  low  of  0.64  to  a  high  of  1.15.  Although  the  variation  was  less  than  the 
zero-crossing  values,  the  values  were  lower  than  desired  for  the  longer  period 
cases  at  0  deg. 

115.  The  shape  of  the  measured  frequency  spectra  showed  excellent 

agreement  with  the  predicted  values  for  the  1.0-sec  cases  regardless  of  wave 
height,  direction,  or  spectral  peakedne.ss.  The  measured  shapes  for  the  longer 
period  cases  did  not  agree  as  well.  The  spectral  pc  iks  were  usually  smaller 
than  the  predicted  values  by  a  factor  of  the  order  of  2.  This  explains  the 
low  values  for  the  heights.  Low  and  high  frequency  tails  were  reason¬ 

ably  well  predicted.  The  RAO  transfer  functions  will  be  used  prior  to  future 
runs  to  compensate  for  these  observed  differences  in  period,  height,  and 
shape . 

116.  The  calculated  peak  periods  and  .  .ean  wave  directions  from  the 
unidirectional  spectral  analysis  agreed  very  well  with  target  values. 

Excellent  agreement  was  reached  between  the  measured  and  predicted  directional 
spreading  estimates  for  the  1.0-sec  cases.  The  longer  period  cases  indicated 
a  wider  spread  than  predicted  which  Increases  with  increasing  period.  Part  of 
this  departure  may  be  explained  by  a  frequency  dependence  of  the  simulation 


and  analysis  software  due  to  random  number  seeds,  truncation  of  the  fourier 
series,  and  sparsity  of  gages.  Analysis  of  these  cases  with  an  MLM  or  MEM 
should  give  higher  resolution  and  narrower  spreading  function  estimates. 
Another  source  of  error  could  be  the  smaller  than  desired  frequency  spectra  at 
the  peak  period.  All  potential  sources  of  error  should  be  investigated 
further . 

117.  Wave  group  statistics  calculated  include  correlation  coefficients, 
probability  distributions,  mean  and  total  run  lengths,  and  correlations  with 
Goda's  peakedness  parameter.  The  calculated  correlation  coefficients  Rhh(k) 
generally  agreed  with  previously  observed  values,  except  for  the  low  1.5-  and 
2.0-sec  cases  at  0  deg.  As  expected,  correlation  drops  off  rapidly  for  lags 
greater  than  one  and  increases  as  the  spectral  peakedness  (i.e.  narrowness)  of 
the  specur^  increases.  The  30-deg  cases  had  uniformly  higher  values  than  the 
corresponding  0-deg  cases,  indicating  that  a  trigonometric  relationship  prob¬ 
ably  exists.  A  correlation  for  every  fifth  wave  height  was  indicated  for  the 
1.5-  and  2.0-sec  cases.  These  trends  should  be  investigated  further  if 
possible . 

118.  For  the  Hm  threshold,  the  Kimura  model  better  matches  the 
observed  probabilities.  For  the  H^g  threshold,  the  match  is  not  as  good  but 
still  better  than  that  for  the  Goda  model.  As  expected,  the  Goda  model  better 
predicted  the  broader  spectral  cases  with  Rhhd)  1  ess  than  0.2. 

119.  Mean  run  lengths  j1  predicted  by  the  Kimura  method  are  slightly 
higher  than  measured  values.  Goda's  model  predicts  a  constant  value  which 
seriously  underp^ed icts  all  cases  except  for  broad  spectral  cases  with 
Rhh(1)  >  0.2  .  The  30-deg  cases  had  slightly  higher  values  than  the  corre¬ 
sponding  0-deg  cases. 

120.  Total  run  lengths  show  a  strong  trend  for  both  thresholds  for  both 
0-  and  30-deg  cases.  The  1^  threshold  trends  show  less  scatter  than  the 
Hmg  threshold  values.  A  least  squares  fit  could  easily  be  calculated  for 
both  thresholds. 

121.  Although  Goda's  spectral  peakedness  parameter  Qp  is  probably  not 
the  best  indicator  of  wave  groupiness,  it  does  show  general  trends.  Increases 
in  Qp  are  observed  with  increases  in  J,  and  R^(1)  •  strong  trend  is 
evident  in  the  0-deg  cases  for  both  thresholds.  For  the  30-deg  cases,  more 
scatter  is  present  for  Qp  >  9  .  Finally,  R^(  1 )  appears  to  level  off  for 
Qp  >  5.5  . 
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APPENDIX  A: 
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m 

Slope  term  of  0 

m 

Overall  mean  wave  direction  for  all  frequencies 
Correlation  parameter 
3 . 14159 

Directional  spreading  standard  deviation 
Left  spectral  width  parameter 
Right  spectral  width  parameter 
Mean  spreading  standard  deviation 
Constant  term  of  a 


a.  Slope  term  of  a 

1  K  m 

4>  Directional  spectra  dummy  variable 

9  Finite  depth  spectral  factor 

Independent,  random  phase 

TMA  finite-depth  parameter 


